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Abstract 

The Pauli operator describes the energy of a nonrelativistic quantum particle with 
spin i in a magnetic field and an external potential. A new Lieb-Thirring type inequality 
on the sum of the negative eigenvalues is presented. The main feature compared to earlier 
results is that in the large field regime the present estimate grows with the optimal (first) 
power of the strength of the magnetic field. As a byproduct of the method, we also obtain 
an optimal upper bound on the pointwise density of zero energy eigenfunctions of the 
Dirac operator. The main technical tools are: 

(i) a new localization scheme for the square of the resolvent of a general class of 
second order elliptic operators; 

(ii) a geometric construction of a Dirac operator with a constant magnetic field that 
approximates the original Dirac operator in a tubular neighborhood of a fixed field line. 
The errors may depend on the regularity of the magnetic field but they are uniform in 
the field strength. 
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1 Introduction 

1.1 Notations 

Let B e C 4 (R 3 ; R 3 ) be a magnetic field, div B = 0, and V € L^R 3 ) a real valued potential 
function. Let A : R 3 — > R 3 be a vector potential generating the magnetic field, i.e. B = Vx A. 
The 3-dimensional Pauli operator is the following operator acting on the space of L 2 (R 3 ; C 2 ) 
of spinor- valued functions: 

H = H{h, A, V) := [tr ■ (-ihV + A)] 2 + V = (-ihV + A) 2 + V(x) + ha ■ B(x) , (1.1) 

where a = (a 1 , a 2 , cr 3 ) is the vector of the Pauli spin matrices, i.e., 




The spectral properties of H depend only on B and V and do not depend on the specific 
choice of A. We shall be concerned only with gauge invariant quantities therefore we can 
always make the Poincare gauge choice. In particular, we can always assume that A is at 
least as regular as B. The operator H = H(h, A, V) is defined as the Friedrichs' extension of 
the corresponding quadratic form from C^°(R 3 ; C 2 ). 
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The Pauli operator describes the motion of a non-relativistic electron, where the electron 
spin is important because of its interaction with the magnetic field. For simplicity we have 
not included any physical parameters (i.e., the electron mass, the electron charge, the speed 
of light, or Planck's constant h) in the expressions for the operators. In place of Planck's 
constant we have the semiclassical parameter h and in most of the paper we also set h — 1. 

The last identity in can easily be checked. If we define the three dimensional Dirac 
operator 

V := a ■ (-ihV + AO)) , (1.2) 

then we recognize the last identity in as the Lichnerowicz' formula. 

The eigenvalues of H below the essential spectrum are of special interest. They determine 
the possible bound states of a non-relativistic electron subject to the magnetic field B and 
the external potential V. Under very general conditions on V and B one can show that the 
bottom of the essential spectrum for the Pauli operator is at zero (see [HNW ). This is in 
sharp contrast to the case of the spinless magnetic Schrodinger operator, (— ihV + A) 2 + V(x), 
whose essential spectrum is not known in general even for decaying potentials. 

Therefore we shall restrict our attention to the negative eigenvalues, &i{H) < 62(H) < 
. . . < of H. It is known that under very general conditions there are infinitely many 
negative eigenvalues even for constant magnetic field |Solj . [Sob-86j, however their sum is 
typically finite. We recall that the sum of the eigenvalues below the essential spectrum is 
equal to the ground state energy of the noninteracting fermionic gas subject to H. 

The sum of the negative eigenvalues, ^ . ej (H) , has been extensively studied recently. 
In order to find the asymptotic behavior of the ground state energy of a large atom with 
interacting electrons, one needs, among other things, a semiclassical asymptotics for ^ . ej(H) 
as h — > 0. 

The semiclassical formula for the sum of the negative eigenvalues is given as 

E scl (h,B,V) := -h- 3 [ P(h\B(x)\,[V(x)]-)dx (1.3) 

JR3 



with 

B ( 00 \ 9 00 

P(B, W):= — \ W 3/2 + 2 V [2uB - Wf! 2 ) = — V d v B[2vB - Wf! 2 (1.4) 

being the pressure of the three dimensional Landau gas (B, W > 0). Here [x]_ = max{0, — x] 
refers to the negative part of x, d := {2ir)~ 1 and d u := tt~ 1 if v > 1. Observe that if = 
o(/i _1 ) then E sd reduces to leading order to the standard Weyl term, — 2(157r 2 )~ 1 /?r 3 / R 3|V]_ , 
as h — > 0. The main feature of the semiclassical formula is that it behaves linearly with the 
field strength in the strong field regime. 
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For the proof that X]j e j(^0 * s asymptotically equal to E sc i as h — * 0, first one must 
establish a non-asymptotic bound on the sum of the negative eigenvalues to control various 
error terms from the non-semiclassical regions. Such estimates for general Schrodinger type 
operators are often referred to as Lieb-Thirring (LT) type estimates fLTlJ. The bound must 
behave like the semiclassical formula in all relevant physical parameters; in this case, in 
particular, it should grow linearly in the field strength. A weaker apriori estimate typically 
leads to a semiclassical asymptotics that is not uniform in the field strength jSob-98j , jES-IIJ . 



1.2 Summary of previous results 

A non-asymptotic LT bound for the Pauli operator has first been established in jLSY-Hj for 
the case of the constant magnetic field, B = const., 

\ej(H)\ < (const.)( j [vf 2 + f |B|[F] 3 _ /2 ) (1.5) 

3 J J 

with h = 1 and this bound was used to prove that E sc i gives the correct asymptotics for the 
sum of the negative eigenvalues. 

The first generalizations of such estimates for non-homogeneous magnetic fields were given 

in |E-1995j . The first general bound was of the form (const.) ( J[V] 5 J 2 + ||B||^ 2 /[V]_ j , then 

the main focus was to study unbounded fields. It was observed, that (jl.5|) cannot hold in 
general. There are two problems in connection with (jl.5j) for nonhomogeneous field. 

Firstly, even when B has constant direction in R 3 (J1.5)) is correct only if |B(x)| is replaced 
by an effective field strength, B e ^(x), obtained by averaging |B| locally on the magnetic 
lengthscale, |B|~ 1//2 . 

Secondly, the existence of the celebrated Loss-Yau zero modes |LYj contradicts (jl.5[) . In- 
deed, for certain magnetic fields with nonconstant direction the Dirac operator T> has a non- 
trivial L 2 -kernel. In this case a small potential perturbation of T> 2 shows that |e,, (if)| 
behaves as j n(x)[V(x)]_dx, i.e. it is linear in V. Here n(x) is the density of zero modes, 
n(x) = \ u j{ x )\ 2 ) where {uj} is an orthonormal basis in KerP. Thus an extra term linear 
in V must be added to (jl.5|) . It turns out that in order to estimate n(x) by the magnetic field 
it is again important to replace |B(x)| by an effective field. 

The problem of the effective field was first succesfully addressed by Sobolev, [Sob-96j, 
Sob-97j and later by Bugliaro et. al. BFFGS and Shen [HE]- In particular, the L 2 -norm 



of the effective field, ||-B e ff||2, i s comparable to ||B|| 2 in |BFFGSj . and the same holds for any 
L p -norm in Shen's work. In a very general bound proved in |LLSj the second term in (jl.5|) is 
replaced with ||B||2^ 2 || V^. 
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In the works IE-1995] . |Sob-9 7 . [Shj, |LLSj . jBFFGSj on three dimensional magnetic Lieb- 
Thirring inequalities, the density n(x) is estimated by a function that behaves quantitatively 
as |B(x)| 3 / 2 . In particular, in the strong field regime these estimates are not sufficient to 
prove semiclassical asymptotics uniformly in the field strength, they typically give results up 
to HBlloo < (const.)h- 1 |Sob-98j . 

We remark that the bounds in |LLSj and [BFFGS have nevertheless been very useful in 
the proof of magnetic stability of matter. In this case the magnetic energy, J |B| 2 , is also part 
of the total energy to be minimized, therefore even the second moment of the magnetic field 
is controlled. We also remark that if the field has a constant direction, then no Loss-Yau zero 
modes exist, n(x) = 0. In this case Lieb-Thirring type bounds that grow linearly with |B| 
have been proved in |E-1995j and [Sob-96] . |Sob-97j . This problem is technically very similar 
to the two dimensional case. 

Since n(x) scales like (length) -3 and |B(x)| scales like (length)' 2 , a simple dimension 
counting shows that n(x) cannot be estimated in general by the first power of |B(a;)| or by 
any smoothed version B c g(x). However, if an extra lengthscale is introduced, for example 
certain derivatives of the field are allowed in the estimate, then it is possible to give a bound 
on the eigenvalue sum that grows slower than |B| 3 / 2 in the large field regime. There are only 
two results so far in this direction. 

The work |BFGj uses a lengthscale on which B changes. The estimate eventually scales 
like 6 17//12 , if the magnetic field is rescaled as B(x) i— > 6B(x), b 3> 1. As far as local regularity 
is concerned, only B G H} oc is required. However, n(x) is estimated by a quantity that 
depends globally on B(x) not just in a neighborhood of x. On physical grounds one expects 
the following locality property: the zero modes of T> are supported near the support of the 
magnetic field. 

We prove a stronger locality property, namely that the size of |B| away from a compactly 
supported negative potential will be irrelevant for the estimate on the sum of the negative 
eigenvalues. The result of |BFGj does not give such bound for an important technical reason. 
In order to produce an effective field strength B c r, the |B| is averaged out by a convolution 
function ip that must satisfy |V</?| < (const.)ip, i.e. tp must have a long tail. For B e L 2 
the effective magnetic field has a comparable L 2 -norm, but it is not true for the localized 
L 2 -norms. 

Our earlier work |ES-lj had a different approach to reduce the power 3/2 of |B| in the 
estimate of n(x). We introduced two global lengthscales, L and t respectively, to measure the 
variation scale of the field strength |B| and the unit vector n := B/|B| that determines the 
geometry of the field lines. This required somewhat more regularity on B than [BFG] and it 
also involved the unnatural H /1,1 -norm of V. The estimate behaved like 6 5//4 in the large field 
regime, if we rescaled B i— ► 6B, b ^> 1. For fields with a nearly constant direction, I ^> 1, 
the bound was actually better, it behaved like b + 6 5//4 £~ 1//2 . This indicates that it is only the 
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variation of n and not that of B that is responsible for the higher 6-power. 

Due to the improvement in the 6-power from 3/2 to 5/4 in the Lieb-Thirring estimate we 
could also prove the semiclassical eigenvalue asymptotics in the regime b <C h~ 3 for potentials 
in W 1 ' 1 jES-IIj . This bound turned out to be sufficient to show that the Magnetic Thomas- 
Fermi theory exactly reproduces the ground state energy of a large atom with nuclear charge 
Z in the semiclassical regime, i.e. where b <C Z 3 , Z — > oo [ES-IIj. The condition b <C Z 3 
is optimal as far as the semiclassical theory is applicable as the results of [LSY-I show for 
super-strong (b > Z 3 ) constant magnetic fields. 

Despite the successful application of the bound in [BFGj to the stability of matter with 
quantized electromagnetic field with an ultraviolet cutoff [BFrO] . and despite that the Lieb- 
Thirring inequality given in |ES-lj fully covered the semiclassical regime of the large atoms, 
it is still important to establish a uniform Lieb-Thirring type bound with the correct power 
in the magnetic field. Such bound will likely be the key to generalize the analysis of the 
super-strong field regime of |LSY-Ij to non-homogeneous magnetic fields. In this paper we 



present a Lieb-Thirring bound that 

• grows linearly in the field strength; 

• depends on the potential V in a natural way; 

• has the locality property in the sense discussed above. 

We also state the corresponding semiclassical result in Theorem 13.31 but its details, that are 
similar to |ES-IIj . will be published separately. 

A simpler proof of a Lieb-Thirring estimate with both the linear dependence in the field 
strength and the correct behavior in V is given in jES-IVj . This approach, however, does not 
give the locality property. 



1.3 Density of zero modes 

As a byproduct, we also obtain a bound on the density of the zero modes, n(x), that behaves 
optimally in the field strength in case of regular fields. Actually, we control the density of 
all low lying states by giving an estimate for the diagonal element of the spectral projection 
kernel U(V 2 < c)(x,x) that grows linearly with the strength of the magnetic field for any 
fixed constant c. 

We remark that the zero modes of the Dirac operator for particular classes of magnetic 
fields are well understood. The surprising first examples were due to Loss-Yau |LY| and later 
the present authors gave a more systematic geometric construction [ES-III . This construction, 
in particular, gives examples that show that the density can grow at least linearly in the 
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magnetic field strength. Other generalizations of the original construction of Loss-Yau are 
also available [AMNj . |El-lj . However, there is no complete understanding of all magnetic 
fields with zero modes yet. 

It is also known that magnetic fields with zero modes form a slim set in the space of all 
magnetic fields ( jBEj . |El-2j ) but no quantitative result is available in the general case. Our 
result (Corollary I3.2|) is the first general estimate on the density of zero modes that scales 
optimally (i.e., linearly) in the field strength. This result is formulated as a corollary, since it 
easily follows from the main theorem, but we shall prove it first on the way to the proof of 
our main theorem (Theorem 13. 1|) . 

It is amusing to note that it takes a considerable effort to show that zero modes exist at 
all, but it is even more difficult to give an optimal upper bound on their densities for strong 
regular fields. This is actually the main technical achievement of the present paper. 

1.4 Organization of the proof 

In Section |21 we introduce what we call the combined lengthscale L c (x) of a given magnetic 
field B(x). This is a local variation lengthscale on which the magnetic field does not change 
substantially. More precisely, this is the case only in the regime where the magnetic field is 
strong; where the field is weak, the combined lengthscale is simply chosen to be of the order 
of the magnetic lengthscale, |B| -1 / 2 . 

In Section |3] we formulate our main result on the new Lieb-Thirring inequality (Theorem 
13. 1|) and its corollary on the density of zero modes. We also state a semiclassical result 
( Theorem 13. 3|) whose proof will be published separately. 

The proof starts in Section 0] with a separation of the contributions from the low and the 
high energy regimes. The cutoff threshold is space dependent, it is at a level P(x) ~ L c (x)~ 2 . 
Technically it is done by inserting P(x) into the resolvent in the Birman-Schwinger kernel 
and using a resolvent expansion. We will call the two regimes the zero mode regime and 
the positive energy regime, respectively, because the separation is dictated by the need for 
a special treatment of the zero modes. The basic estimates on the contribution from these 
regimes are given in Theorem 14.31 We remark that to ensure ultraviolet convergence in the 
zero mode regime, squares of resolvents need to be estimated as well ( |BFFGS] ). 

In both regimes we perform a two-scale localization, like in |ES-lj . For both localizations, 
however, the approaches used here are substantially improved, as we explain below. 

The first localization is isotropic and its lengthscale is determined by L c (x). This is 
constructed in Section The main difference between the current isotropic localization and 
the corresponding one in ES-Ij is that in our earlier paper we assumed a universal positive 
bound on the combined lengthscale, therefore we could use a regular grid of congruent cubes. 
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In order to ensure the locality property, in this paper we need to use a covering argument to 
select localization domains of different sizes and with a finite overlap. In domains where the 
magnetic field is relatively weak (|B| < (const.) L~ 2 ), we shall neglect all magnetic effects. 

In Sections El and [7| we show how to localize the eigenvalue estimates onto the isotropic 
domains. In the positive energy regime we apply a version of the IMS localization formula for 
the resolvent ( Proposition 16 .lj) that was already used in [BFFGS . However, the same formula 
does not hold for the square of the resolvent which is needed in the zero mode regime. A 
new localization scheme is developed in Proposition 17. II to localize the square of the resolvent 
of a second order elliptic operator. The localized versions of the necessary estimates in the 
positive energy regime and in the zero mode regime are stated in Propositions 16.21 and 17.21 
respectively. 

Typically it is not hard to localize resolvents of second order elliptic operators onto cubes 
of size I at the expense of an error i~ 2 . However localizing the square (or higher powers) 
of the resolvent requires off-diagonal estimates on the resolvent kernel (see Proposition 7.1). 
While these are typically easily available for scalar elliptic operators without spin, we do not 
know any apriori off-diagonal control on the resolvent of D 2 . If the original Pauli operator 
is estimated by a constant field Pauli operator, then aposteriori we can extract off-diagonal 
estimates, but without comparison with the constant field problem, we do not have off-diagonal 
control. This is the main reason why we are unable to extend the elegant and short method 
of jES-IVj to give any locality properties. 

Starting from Section |H] a second localization is performed onto curvilinear cylindrical 
domains with a transversal lengthscale |B(x)| -1 / 2 along the field lines. The geometry of 
the cylindrical domains and the coordinate system are explained in Section 18.11 and a new 
partition of unity subordinated to the cylindrical domains is constructed in Section IHT^l 

Within each cylindrical domain the magnetic field is approximated by a field /3 C , given as a 
2-form, that is constant in the appropriate cylindrical coordinates and after a conformal change 
of the metric (Definition 19. 1|) . The Dirac operator T> c with a magnetic field /3 C (Definition ^. 3j) 
will be used to approximate the original Dirac operator T> in the corresponding cylindrical 
domain. Section El is devoted to the construction of V c and it uses the geometric structure 
behind the Dirac operator on a non-flat manifold outlined, for example, in jES-IIIj . 



The second main difference between [ES-Ij and the current work lies in the cylindrical 
localization. In [ES-IJ we considered straight cylinders to approximate tubular neighborhoods 
of magnetic field lines and we approximated the field by a constant one within each cylinder. 
The curving of the magnetic field was not respected by the approximation hence the error was 
not uniform in the field strength. This is the main reason why the Lieb-Thirring inequality 
in |ES-Ij does not have the optimal |B|-power. 

In the new construction the cylindrical localization domains are curved in such a way as 
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to follow a field line and we also construct appropriate spinor coordinates. This geometric 
approach enables us to control T> 2 with errors that are uniform in the field strength although 
they depend on the combined lengthscale of B. This eliminates the |B|-dependent error in 
the large field regime. The near-zero energy states, in particular the zero modes of T>, need 
to be controlled with such a precision in order not to overestimate their contribution to the 
negative eigenvalues of T> 2 — V . The proof of the |B|-independent control is quite involved 
and it relies heavily on the intrinsic geometric properties of the Dirac operator. 

Section El completes the proof of the positive energy regime. In this regime errors that 
are independent of |B| can be absorbed into the local energy shift P(x). Proposition 110. ll 
contains the necessary spectral estimate localized onto the cylindrical domains in the original 
coordinates. We first translate the estimate into cylindrical coordinates. In these coordinates 
the approximating field is constant and we can use the magnetic localization formula (|10.18|) 
from |ES-IIj . This method yields a |B|-independent cylindrical localization error. Finally, 
having constructed the approximating local Dirac operators with constant fields, we can use 
the Lieb-Thirring inequality for the Pauli operator with a constant field obtained in |LSY-II| . 

In Section ^2 we complete the estimate of the zero mode regime. We need to estimate 
the density of the near-zero energy states of T> 2 . This is given by the diagonal kernel of the 
spectral projection operator, Ii{V 2 < P )(x,x) where P is the typical value of the regularly 
varying function P around x. This operator can be bounded by the resolvent, but the diagonal 
element of the resolvent is infinite because of the ultraviolet divergence. Therefore we need 
to control II {T> 2 < P ) by the square of the resolvent, [T> 2 + Pq)~ 2 . F° r regions with weak 
magnetic fields the magnetic field can be neglected and we can simply use the diamagnetic 
inequality. The problem thus can be reduced to estimating the resolvent square of the free 
Laplacian ( Section II 1.1)1 . 

For regions with a strong field (Section lll.2)) we again use the approximating constant field 
operators. However, the magnetic localization formula ()10.18j) is not valid for the square of 
the Pauli operator, so localizing onto cylindrical domains is more complicated. Fortunately, 
at this stage we do not need operator inequalities, we need to estimate only the diagonal 
element of the square of the resolvent at each fixed point u. First we transform the problem 
into the new coordinates associated with the field line through u (estimates (jll.6)) and (JTTTTJ)). 
Then we use resolvent expansions extensively to approximate T> by T> c . Since we estimate the 
square of the resolvent, we need to control the off diagonal elements of the resolvent itself. 
For a constant field, the offdiagonal decay is Gaussian on a magnetic lengthscale |B| -1 / 2 . 
A similar feature is proved for the resolvent with a nonconstant field via the constant field 
approximation. 

In the Appendix we collected the proofs of several Propositions and Lemmas which can 
be skipped at a first reading. 

It is amusing to note that the most complicated part of the proof (Sections [7| and fTT|) 
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controls the possible ultraviolet regime of near zero energy states. On physical grounds this 
regime should be irrelevant if we knew that low energy eigenstates of T> 2 have transversal 
momentum of order |B| X / 2 and parallel momentum independent of the field strength. The 
main difficulty is to obtain such information on the low lying states. 

Convention: Throughout the proof universal constants are denoted by a general c whose 
value can be different even within the same equation. Constants depending on numbers a,b, . . . 
are denoted by c(a, b, . . .). Integration over R 3 with respect to the Lebesgue measure, f R3 dx, 
is simply denoted by J. We shall say that two positive numbers a, b are comparable if 
\ < a/b < 2. 



2 Lengthscales of the magnetic field 

Let B e C 4 (R 3 ,R 3 ) be a magnetic field and let n := B/|B| be the unit vectorfield in the 
direction of the magnetic field at all points where B does not vanish. 
For any L > and x G R 3 we define 

B L (x) :=sup{|B(y)| : \x - y\ < L) (2.1) 

and 

b L (x) :=inf{|B(y)| : \x — y\ < L} (2.2) 

to be the supremum and the infimum of the magnetic field strength on the ball of radius L 
about x. These functions are continuous in both the L and x variables. 

The Pauli operator will be localized on different lengthscales determined by the magnetic 
field. We now define these scales. 

Definition 2.1 (Lengthscales of a magnetic field). Given a C 4 -magnetic field B. We 
define the magnetic lengthscale of B as 

L m (x) := sup{L > : B L (x) < L~ 2 } . (2.3) 

The variation lengthscale of B is given by 

L v (x) := supjL > : L 7 sup 1 1 V 7 B(y) : \x - y\ < l| < b L (x), 7 = 1, 2, 3, 4| (2.4) 

Finally we set 

L c (x) := m&x{L m (x) , L v (x)} (2.5) 
to be the combined lengthscale o/B at x. 
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A magnetic field B : R 3 — > R 3 determines two local lengthscales. The magnetic length- 
scale, L m , is comparable with |B| -1 / 2 . The lengthscale L v determines the scale on which the 
field B varies. One may think of L v as the smaller of the two lenghtscales describing the 
variation of the field strength, i.e., the variation scale of log |B|, and the variation scale of the 
field lines n. 

For weak magnetic fields the magnetic effects can be neglected in our final eigenvalue esti- 
mate, so the variational lengthscale becomes irrelevant. This idea is reflected in the definition 
of L c ; we will not need to localize on scales shorter than the magnetic scale L m . 

Note that for any B G C 4 (R 3 ,R 3 ) we have that < L c (x) < oo for all x G R 3 . If 
L c (x) = oo for some x G R 3 , then B is constant on R 3 . Moreover the value L c (x) at any x 
does not depend on B outside the ball centered at x with radius L c (x). If follows in particular, 
that if B vanishes in a ball of radius 5 around x, then 5 < L c (x). 



3 Main Theorem 

We are ready to state our main results. 

Theorem 3.1 (Uniform Lieb-Thirring inequality). We assume that the magnetic field 
is B G C 4 (R 3 ,R 3 ). Let A G C 4 (R 3 ,R 3 ) be a vector potential, V x A = B, and let V : = 
cr ■ (— i V + A) be the free Dirac operator with magnetic field B on the trivial spinorbundle over 
R 3 ; that can be identified with L 2 (R 3 ,C 2 ). Let V be a scalar potential. Then the sum of the 
negative eigenvalues, Cj, of the Pauli operator H := V 2 + V = [cr ■ (— iV + A)] 2 + V satisfies 

\TrH_\ = J2\e J \<cJ[V] 5 / 2 + c J \B\[VfJ 2 + cJ(\B\ + Lf)L^[V}_ (3.1) 

with universal constants. 

Notation: For any self-adjoint operator H we let H := j[\H\ — H] denote its negative 
part. 

Corollary 3.2 (Density of zero modes). Given a magnetic field B G C 4 (R 3 ,R 3 ) with a 
combined lengthscale L c , the density of zero modes of the free Dirac operator!) with magnetic 
field B satisfies 

n(x) := M*)| a < c(|B(x)| + Lf(x))L;\x) (3.2) 

3 

with a universal constant, where {uj} is an orthonormal basis in the kernel ofV. 
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Remarks, (i) The density function n(x) was also estimated in |BFG]. In the strong field 
limit B i — ► 5B, b ^> 1, the estimate behaved as b 17 ^ 12 . Moreover, unlike in [BFG] , our estimate 
on n(x) uses only local information on B(x) as explained in Section |21 For example, if B 
vanishes inside a ball centered at x with radius 5, then n(x) < c5~ 3 . 

(ii) The bound ()3.2|) is optimal as far as the strength of the field |B| is concerned. This fact 
follows from the construction of Dirac operators with kernels of high multiplicity following the 
method of |ES-lllj . For example, the density of Aharonov-Casher zero modes for a constant 
magnetic field of strength B ^> 1 on S 2 is of order B. The geometric procedure of [ES-IIIj 



allows one to construct a Dirac operator on R 3 whose zero energy eigenf unctions are obtained 
from the eigenf unctions on S 2 by an explicit transformation. The density of these states 
remain comparable to the strength of the magnetic field at least away from infinity. 

(iii) Notice that the Lieb-Thirring inequality of |LSY-llj for a constant field is recovered 
in Theorem 13.11 

(iv) The uniform Lieb-Thirring bound for a constant direction field, [Sob -97j , |ES-Ij , does 
not directly follow from our main theorem as it is stated. On one hand, ()3.1|) contains a term 
linear in V that is unnecessary for a constant direction field. On the other hand, we assume 
high regularity on B. This regularity is needed only to construct the appropriate curvilinear 
cylindrical localization, which is unnecessary for a field with constant direction. 

However, our present technique to estimate squares of the resolvents can improve these 
results in another aspect. For example, if the support of B and the support of V are sepa- 
rated, our Lieb-Thirring estimate depends only on the separation distance whereas all previous 
bounds scale with the magnitude of B. As a byproduct of such a result one can also improve 
the estimates on the ground state density of the two dimensional Pauli operator given in 
jF93j . 

Armed with a uniform Lieb-Thirring inequality, the following semiclassical asymptotics 
may be proved by combining the techniques of the current paper and [ES-II] . The details of 
the proof will be published separately. 

Theorem 3.3. We assume that B G C 4 (R 3 ,R 3 ) and V G L 5/2 (R 3 ) n L^R 3 ). Then the 
sum of the negative eigenvalues, ej(b,h), of the Pauli operator [a ■ (—ihV + bA)] 2 + V is 
asymptotically given by 

l im ^ ; = i (3.3) 

h-o E scl {h, bB, V) 1 ; 

where the limit is uniform in the field strength b. 



Remark. This result was obtained for a homogeneous magnetic field in |LSY-IIj . Analogous 



results for d — 2 were obtained in jES-IIj and |Sob-98j . The latter work also extends the 
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two dimensional analysis to obtain (|3.3|) for three dimensional magnetic fields with constant 
direction. For a general three dimensional magnetic field the limit (|3.3|) is proven up to 
b <C h~ 3 for V G W 1 ' 1 in ES-II . With a different method Sobolev also obtains ()3.3|) up to 
b < (const. )h~ 1 without assumptions on the derivatives of B and V [Sob-98j. 



4 Proof of the Main Theorem 13.11 
4.1 Tempered lengthscale 

Since localization errors decrease with the localization length, we would optimally like to 
choose the biggest possible scale, i.e. L c (x), for our localization scale. However, neighboring 
localization domains must be comparable in size so that the localization errors could be 
reallocated. This forces us to require a tempered behavior on the localization scales, which 
may result in choosing a localization scale smaller than L c . Proposition 14.21 below shows 
that this technical requirement can be met at the expense of a factor | and this justifies the 
introduction of the tempered lengthscale L := \L C . Before the precise statement we need the 
following definition: 

Definition 4.1. Let e > be a positive number. A positive function f(x) on R 3 is called 
^-tempered if 

\x-y\<e~ 1 f{x)^\< S -^\<2 \/x,yeR 3 . (4.1) 
If e = 1, then a 1 -tempered function will be simply called tempered. 

Proposition 4.2 (Existence of tempered lengthscale). For any not identically constant 
magnetic field, B G C 4 (R 3 ,R 3 ) ; L(x) = \L c (x) is finite and defines a tempered function. 
Moreover, if B L ^(x) > L(x)~ 2 , then b L ^(x) > 0, and for 7 = 1,2,3,4 

L(x) 7 sup{|v 7 |B(y)| J : \x-y\<L(x)}<b L(x) (x), (4.2) 

and 

L(x) 7 sup{|V 7 n(?/)| : \x - y\ < L(x)} < 1. (4.3) 
For a constant magnetic field B = const we have L c = 00 and we set L(x) := 00. 

The proof is given in Section IA.1I 

For a constant magnetic field B the tempered scale has been defined to be infinity for the 
transparent formulation of our theorem. The estimate (J3.1)) for this case has been proven in 
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|LSY-llj . It is possible to apply our proof to this case as well, but setting L = oo directly may 
require minor remarks along the proof. In order to avoid this inconvenience, we can choose L 
to be any fixed real number for which the proof goes through without changes and finally let 
L — > oo in the final result (|3.1|) . 

We introduce a universal constant < e < that has to be chosen small enough for the 
proof to work but we shall not keep track of the exact numerical value needed. We consider 
it fixed throughout the proof. 

Let L(x) be the tempered lengthscale of B. Introduce £(x) := sL(x), then the properties 
of L(x) set in Definition 12 .11 and Proposition 14.21 are translated into I as follows: 

• t(x) is ^-tempered. 

• For all x G R 3 such that sup{|B(y)| : \x — y\ < e~ x t} > e 2 £~ 2 and for 7 = 1, . . . , 4 we 
have 

£(x) 7 sup{ V 7 |B(y)| : \x-y\ < e -1 *(a;)} < e 7 inf {|B(y)| : \x - y\ < e -1 ^(a:)} . 

(4.4) 

and 

£(x) 7 sup| V 7 n(?/) : \x-y\ < < e 7 . (4.5) 

We define 

P(x) := e- 5 £(x)- 2 
and for any positive function / > we introduce the notation 

R, = R(f) = (V 2 + f)- 1 . 
4.2 Separation into low and high energy regimes 

We shall prove Theorem 13.11 with L c replaced by L. Since T> 2 + V > T> 2 — [V]_, we may 
consider only non-positive potentials. For convenience we change the sign and we will work 
with the operator H := T> 2 — V with V > 0. 
By the Birman-Schwinger principle 



\TtH_\ = ^ n(v 1/2 R E V 1/2 ,l^jdE 



(4.6) 



where n(A, fx) is the number of eigenvalues of the operator A greater than or equal to fi. For 
any E > we have, by the resolvent identity, that 

Re — Rp+E + Rp+ePRe — Rp+E + Rp+ePRp+e + Rp+ePRePRp+e ■ 
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Using that P < V 2 + P + E and R E < E' 1 , we obtain 

Re < 2Rp + e + E 1 Rp + eP 2 Rp+e ■ 
For any positive operators Xi,X 2 , 

n{X 1 + X 2 , et + e 2 ) < n{X u e x ) + n{X 2 , e 2 ) , (4.7) 
hence (|4.6|) is estimated as 

|Tr H-\ < J™ n (v 1/2 R P+E V 1/2 , \\ dE + J™ n ^2V 1/2 R p+ eP 2 Rp +e V 1/2 , E^j dE . (4.8) 

The second term carries the contribution of the near zero energy eigenfunctions of the free 
Pauli operator T> 2 . This will be called the zero mode regime. 
For the first term we notice that 

n (v 1/2 R P+E V 1/2 , | J dE = |Tr(P 2 + P - 4V% | (4.9) 

by the Birman-Schwinger principle. This term contains the contribution from free eigenfunc- 
tions with energy at least O(P) and it will be called the positive energy regime. 

The following Theorem estimates the two terms in ()4.8|) and it completes the proof of the 
Main Theorem by choosing e sufficiently small. □ 

Theorem 4.3. For a sufficiently small universal e and with the notations above we have 

Ti(V 2 + P - AV)- < c(e) / (V 5/2 + |B|y 3/2 ) , (4.10) 
jT n^2V 1/2 R P+E P 2 Rp + EV 1/2 ,E^dE < c{e) J VP 1/2 (\B\ + P) . (4.11) 
The proof of Theorem 14.31 is given in the rest of the paper. 

Convention about operator kernels: If A is a Hilbert-Schmidt operator on a Hilbert space 
of the form L 2 (d/z) ® C^, iV e N, we denote by A(x, y) its x A^-matrix valued integral kernel 
which is L? on the product space. If, in addition, A is of trace class, we can even define its 
diagonal kernel which, by a slight abuse of notation, will be denoted by A(x, x). One possible 
way to define it is to write A as a product of two Hilbert-Schmidt operators, A = HK, and 
A(x, x) := J H(x, y)K(y, x)dfi(y). This is an L 1 matrix valued function of x and as such it is 
independent of the choice of H and K. 

Convention about traces: We shall denote by Tr the trace on L 2 (d/x) <S> C N and by tr the 
trace on C N . If A is of trace class on L 2 (dfi) £g> C^, then trA(x,x) is in L 1 (d y u). 
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5 Isotropic geometry of the first localization 



In this section we construct the domains for the first localization. The construction is deter- 
mined by the function £{x). We shall construct a discrete set of points {xi}. The localization 
domains will be balls about xi with radii £{xA and they will have finite overlap. Moreover, the 
magnetic field will not change much in each localization ball since £(x) determines the local 
scale of variation of B. Outside of this domain the field will be replaced by a constant field. 
This procedure will apply to balls with relatively strong fields. On balls where B is small we 
neglect magnetic effects and replace the Pauli operator by the free Laplacian. 



5.1 Regular fields 

Definition 5.1. Given £, K > and a ball D of radius I centered at zq G R 3 . A magnetic 

field B is called .D-strong if |B(^o)| > e~ 2 £~ 2 , otherwise it is called D-weak. A D-strong 
magnetic field is called (D, i^)-regular if for 7 = 1, . . . , 4 



(i) 



V 7 IB| 



< Ke^£-^\B(z )\ on D; 
(11) |V 7 n| < Ke^£-^ on D (with n := B/|B|j; 

A (D,K) -regular field B is called extended (D, i^)-regular if it is continuous on the 
whole space and B is constant outside of D. The value o/B outside of D is denoted by B^ 
and for ^ we set n M := Boo/IB^I . 



A (D, i^)-regular field B clearly has a small total variation on D: 



\B(x) 



(5.r 



|B(z )| <2Ke\B(z )\ 

for any x G D. For an extended (D, i^)-regular field (|5.1|) is valid for any x G R 3 , and 

||n(:r) - n^ll < Ke . (5.2) 
The following statement follows from the definitions above: 

Lemma 5.2. Let B{x) and £(x) satisfy the conditions \4-l\j , \4-4\) an d Let D be the 

ball of radius 10£(zq) about some z G R 3 . 

(i) IfB is D-strong, then B is (D, l)-regular and |B(sc)| > e~ 1 £(zq)~ 2 for any x G D. 



(ii) IfB is D-weak, then |B(x)| < e 



-2/ 



for any x G D. 



Proof, (i) For e < we obtain that for any x G D 

\B(zq)\ - \B(x)\ < lOeinf {|B(y)| : \x - y\ < £ _1 ^ )} < lOe\B(z ) 
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using (J4.4j) . In particular, |B(x)| > (1 — 10e)|B(z )| > e~ 1 £(xy 2 for any x E D because B 
is D-strong, |B(^ ) I — £~ 2 {^£{zo))~ 2 , and £(z$) < 2£(x). Properties (i) and (ii) in Definition 
O follow from (jOj) and (|Ojl . 

(ii) Suppose that for some i 6 Dwe have |B(rc) | > 6~ 2 £(zq)~ 2 . Using that \zo — x\ < 
10£(^ ) < 20£(x) < e-^(a;) if e < ^, we obtain 

|B(^ )| - |B(a;)| < lOeinf {|B(j/)| : |y - ac| < e -1 ^)} < 10e|B(zb)| 

which contradicts to |B(z )| < £~ 2 (l0£(z ))~ 2 . □ 



5.2 Covering lemma and cutoff functions 

Let B(x,r) denote the closed ball centered at x with radius r. We introduce the following 
notations for any x E R 3 

D x :=B(x,^j, D x :=B(xJ(x)^, D x := b(x,1M{x)) . 

Definition 5.3. Let £(x) be an e-tempered function and let I be a countable index set. The 
discrete set of points {xi}i £ j is called an ^-uniform set of points with intersection constant 
N if 

fzjR 3 cU ie A; 

(ii) Any ball D x intersects no more than N other balls from the collection {D x .}. 

The proof of the following covering Lemma is given in Section [PI 

Lemma 5.4. Let £(x) be e-tempered, then there exists an £-uniform set of points {xj} i6 / with 
some universal intersection constant N . 

In the rest of the proof we fix such a collection of points {xi}, determined by the magnetic 
field via £(x). For brevity we shall use £{ := £(xi), Di := D Xi , Di := D Xi and Di := D Xi . 

Definition 5.5. An index i E I, the corresponding point Xi and ball Di are called strong 
(weak) if Ii is Di-strong (weak). 

The following Lemma is an application of Lemma 15.21 

Lemma 5.6. Let Xi be a strong point, then B is (Di, l)-regular and infg |B| > e~ l t[ 2 . If Xi 
is a weak point, then sup^. |B| < e~ 2 £j 2 . 



18 



Given an e-tempered function £(x) and an ^-uniform set of points {xj}j g /, for each i G I we 
choose smooth functions 0j, Xfa Xi an d Xi with values between and 1, such that the following 
hold: 

• E^fO) = 1, supp(0i) C A and ||V0<||oo < < X ; 

• £i = 1 on 5(^,34), supp(£i) C B(x i ,4£ i ), ||Vx*||oo < ^f 1 ; 

• Xi = 1 on B(x i ,4:£ i ), supp(xi) C £(^,54), || Vxdloo < 2£j l ; 

• Xi = 1 on flfo, 64), supp(xi) C 7^), HV^ilU < (24)" 7 , 7 = 1, ... ,4. 

Such choice is possible since the balls Di cover. Notice that V£i is supported on the annulus 

At :^B{x i ,U i )\B{x i ,3£ i ) . 

Finally we choose functions {y?i}ie/ such that <pi = 1 on A*, supp(<£>j) C B{x^ 54) \ 24) 
and |V^| < 2^. 

5.3 Approximate magnetic fields and Pauli operators 

We define approximate vector potentials Aj and magnetic fields Bj := V x Aj, i = 1,2, . . ., 
subordinated to the balls Di. The definition is different for weak and strong indices i G I. 

If % G I is a weak index, then let Aj be the Poincare gauge of B with base point Xj, in 
particular V x A; = B on R 3 and | Aj| < c£i sup^ |B| < ce~ 2 t~ x holds true on Di by Lemma 

15.61 We define Aj := A - (1 - Xi)Aj, then Bj = x«B + Vxi x Aj. Clearly A(x) = Aj(x) for 
all x G B(xi, 64) and 

llBiHoo < csup |B| < ce~ 2 £- 2 (5.3) 

A 

with supp Bj C Di. 

If i G / is a strong index, then Aj is given by the following lemma. 

Lemma 5.7 (Choice of the local field on strong balls). Assume that B is (Di, 1) -regular, 
then there exists a vector potential Aj such that Aj = A on B(xi, 64) and the magnetic field 
Bj = V x Aj satisfies 

Bj(x)=B(x), xeB(xi,Ui) and B^x) = B(ar<) x G R 3 \ B{x u 74) . (5.4) 
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Moreover, B, is extended (A; 100)-regular, in particular for 7 = 1, . . . , 4 



^||V^|B,| 11^ < 100^|B^)| , (5.5) 
^||VX 11^ < lOOe-y , (5.6) 
\Bi(x) -B( Xi )\ < 100e\B( Xi )\ (5.7) 



for any x G R 3 . 

Armed with these definitions of A,;, we define 

Vi := a ■ (-zV + (5.8) 

to be the approximating Dirac operator associated with A- The operator A coincides with 
T> on B(xi, 6£j) because A = Aj in this domain, in particular 

V X i = V iX i • (5.9) 

Proof of Lemma \5. 7| Since B is (Aj l)-regular, from (|5.1j) we obtain that 

|B(x) -B(xi)| < e|B(xi)|, x G A . (5.10) 

Let A* be the Poincare gauge on A of the magnetic field B — B(xi), then V x A* = B — B(a^) 
and 

|Af (x)| <e£i\B(xi)\ (5.11) 

for any x G A- We then define 

Bj := V x (xiAf) + B(xj) (5.12) 

Easy calculations show that this field is (Aj 100)-regular and (|5.5jl - (|5.7jl hold. 

The gauge A* + |B(zj) A (• — xi) generates B, hence 

A = A* + |B(xi) A (• - x,) + 
with some 0j : R 3 — > R. Since (x«Af ) + |B(x.j) A (■ — acj) generates B i; we define 

A, := (XiAf ) + iBfo) A (■ - Xi ) + . (5.13) 

Then V x A, = B ; . □ 
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6 Positive energy regime: proof of (14.101) in Theorem 

We recall the set {xi} constructed in Section O and let 

i t := £(xi), P t := P( Xi ) = e-Hr\ b t := |B(x,)| . 
We also recall that for any positive function / we denote the resolvents by 

R f := R[f] := (P 2 + Z)- 1 . (6.1) 
Note that in general Rf and R g do not commute. For simplicity we also introduce 

(flf + Z)- 1 . (6.2) 

Proposition 6.1 (Pull-up proposition). Let I be a countable index set and let g^, i G I, 

be a family of nonnegative smooth functions such that < Yliei 9i( x ) < 00 f or ever U x £ R 3 - 
Let Ai, i G / be a family of positive invertible self-adjoint operators on L 2 (R 3 , C 2 ). Then 

(Erf)£-^x7(Erf)*E*x*- (63) 

Proof of Proposition \6. 1\ This proof is basically given in [BFFGSJ, we repeat it here for 
completeness. All positive self-adjoint operators below are interpreted as quadratic forms. 
We start with the operator inequality 

J* J ^—^J*J < J*AJ (6.4) 

for any positive self-adjoint operator A and any operator J. 

We define a map J : L 2 (R 3 , C 2 ) i-> 0. L 2 (R 3 , C 2 ) =: H as J : ip ^ {g^}. We define an 
operator A on Ti as A : {ipi} 1— > {Aiipi}. It is easy to check that 



rAJ = ^g i A i g i on L 2 (R 3 ,C 2 ), 

i 

J* J = Y,i9 2 i and that (A)- 1 = Ar l . Thus 

(E«?) s^(E«f) - J'Jyrj" * J'Cavj - e*£« • D 

The following proposition is the localized version of (|4.10J) for strong balls and its proof is 
given in Section fTUl 
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Proposition 6.2 (Positive energy regime). Let D be a ball of radius i and let K > be 

a positive number. Let B be extended (D, K) -regular, let the function < x < 1 fre supported 
on D. Then for any positive numbers M,fi > there exists a constant e(M, K, fi) such that 
for any e < e(M, K, jj) we have 

Tr(D 2 + /2e- 5 r 2 - MxV)_| < c(M, K, e) J (V 5/2 + \B\V 3/2 ^j . (6.5) 

Armed with these two Propositions, we can finish the estimate (I4.1(J|) in Theorem 14.31 
Using the finite overlap property of D^s ( Lemma 15. 4j) . and that Oi < Xi 1> we see that 

l<H(x) :=$>?(*) < N . 

iei 

Moreover, by the localization estimate, 

/ i^i 2 > f E / i^i 2 - ^ E / i p ^i 2 - * E<^> ( v ^> > 

iG/ is/ i£j 

hence 

iei ie/ 

where l(-) is the characteristic function. Using (|5.9|) we may simply replace T> 2 by T> 2 on the 
support of Xi- If e is sufficiently small, we obtain 

V 2 + P + E>^Y1 + P * + E ^ ( 6 - 6 ) 

iei 

using the finite overlap property and that P is comparable to Pi = e~Hf on Di. The resolvent 
can be estimated by 

using Proposition 16.11 Hence, by the Birman-Schwinger priciple 
Tr(£> 2 + P - 4V% = jT n(v 1/2 .RF +E V 1/2 ,;|)d£ 



< 



iei 

< [°° Jy^V^XiHPi + E]^ 2 ,^]^ . (6.7) 
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Here we used 

n(ABA, e) = n{B 1 ' 2 A 2 B 1 ' 2 , e) (6.8) 

for any nonnegative operator B and arbitrary operator A with the choice A = S _1 and we 
estimated S~ 2 < 1. 

We also use a strengthening of (|4.7jl . If the positive self-adjoint operators A, B are disjointly 
supported, i.e., there exists an orthogonal projection II such that UAH = A and (I — H)B(I — 
II) = B, then 

n(A + B,e) = n(A,e) +n(B,e) . (6.9) 

The proof is trivial. 

In order to use Proposition 16.21 we have to pull the summation out in (|6.7|) . We split this 
sum into a few infinite sums so that each contain disjointly supported terms. Since the balls 
{Di} have uniformly finite overlap with constant N (see (ii) of Lemma l5.4j) . there exists a 
partition of the index set I — I\ U I2 U . . . U ijv+i such that if j, j' G for any 1 < A; < N+ 1, 
then fl = 0. Such a partition can be obtained by a greedy algorithm. We order the 
index set / in some way and we put each index one by one into one of the sets. We always put 
the new index into one of the sets where it has no conflict with the indices already put into 
this set. A new index j is said to be in conflict with a previously placed index i if Di ClDj 7^ 0. 
Since every index can have a conflict with at most iV other indices, each index can be placed 
somewhere at each step of the placement. 

Hence, using ()4.7|) first, then ()6.9|) . we have 

N+l 



^ iei ' ^ k=i iei k ' 

N+l , , 

k=l ^ ie/fc ' 
N+l , , 

= EEH yl/2x4 ^ [Pj+ ^ 1/2 'WvTT) 

k=i iel k V ' 
is/ ^ ^ 

so combining this estimate with ()6.7|) and applying the Birman-Schwinger principle in the 
opposite direction, we obtain 

Tr(£> 2 + P - 4V)_| < \ Tl ( V i + p i~ M X 2 V)- (6-10) 
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with M := 16iV(JV + l). 

We then apply Prop osit ion 16 . 21 for each strong index i for the ball D = D^, radius I = 10£i, 
and the magnetic field Bj that is extended (Di,K = 100)-regular. For small enough e we 
obtain 



Tr(D 2 + P % -M X 2 V)- 



< c(e 



M ) [ (V 5 / 2 + |B,|V 3 / 2 ) 
< c(e,M) j (v 5/2 + \B\V 3/2 ^j for i strong, 



where the last inequality follows from (|5.1jl . For the weak indices i we use T> 2 
Aj) 2 + a ■ Bj and \\<r ■ Bj|| < ce 3 Pi (see (J3T3J) ) and we obtain 



-iV + 



< 



Tr ((-^V + A,) 2 -M X 2 ^). 



< c(M) / V*' 2 X . 



Tr(X? 2 + P^Mx 2 V). 



by the usual Lieb-Thirring inequality for magnetic Schrodinger operators without spin. Sum- 
ming up these estimates we obtain from (j6.1()|) that 



|Tr(£> 2 + P-4V)_| < c{e) 



< c(e) j {y^ 2 + \B\V 3 / 2 



again by the finite overlap property of B. This completes the estimate (J4.1(Jj) . □. 



7 Zero mode regime: proof of (14.111) in Theorem 14.3 

The estimate (j4.11|) essentially involves estimating the square of the resolvent of T> 2 . However, 
the analog of Proposition 16.11 does not hold for the square of the resolvent, i.e. 

( E sf) * ( E ^ 9 i) ( E 9 i) < E 

with $(£) = t~ 2 is noi true in general. Here is a 2 by 2 matrix counterexample with gi = g% — 
2 -i/2. 

Ai= 2) ' ^ 2= (? 2) • 
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Without such an inequality, we have to use a resolvent expansion. In addition to the square 
of the localized resolvent, we need to control offdiagonal terms. Such an estimate is given in 
the following Proposition, which is a general statement about squares of resolvents of second 
order differential operators. The proof is given in Section TA. 31 



Proposition 7.1 (Pull-in proposition). Given an e -tempered function £(x) and a function 
F(x) > satisfying 

1 Fix) , . 

for all \x-y\< e- x l{x). SetP( x) = e 5 £(x) 2 . Let A = A-V + B be a first order differential 
operator acting on @ fc L 2 (R 3 ) ; 1 < k < oo, with smooth coefficients, i.e. A(x) is a vector of 
k x k matrices, B(x) is a k x k matrix, all smoothly depending on x. We assume that 

supp(x)||<co. (7.2) 

x 

LetT = A* A. 

Given an t-uniform set of points {xi} i£ j as in Lemma [3^1 Let Xh^ijfi be chosen as 
in Section \5.2l We assume that for every % there exists a first order differential operator 
Ai = Ai-V + B l on fc L 2 (R 3 ) such that 

A = Ai, B = Bi onsupp(xi), (7.3) 



and let T{ = A*Ai. Then there exists an Eq depending only on Co in \7.2j) such that for any 
e < Eq and \i > we have 



-F 



t +p+ ,- T +P + ^ c ^ e '(^ 1 T^ A ^ A 'nT^ + mTp^ e ' (7 ' 4) 

where Fi := sup{F(x) : x G Di} and Pi = P(xi). 

Remark. If A, B are well-behaved, then T looks like an elliptic constant coefficient differ- 
ential operator on short scales. In this case the estimate localizes the square of the resolvent 
in such a way that the diagonal element of the operator kernels on the right hand side of (|7.4j) 
remain finite. This is clear for the second term on the RHS since the estimate is integrable 
in the ultraviolet regime (behaves like p~ A in the momentum p). The first term behaves only 
as p~ 2 , but the supports of Q{ and (fi are well separated, which makes the diagonal element 
finite. 

The diagonal kernels of the localized operators are estimated in the following Proposition 
whose proof is given in Section El 
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Proposition 7.2 (Zero mode regime). Let D be a ball of radius £ > with center z G R 3 
and K > 0. We assume that either 

(i) || B || oo < ce~ 2 £~ 2 and B is supported on the ten times bigger ball D = B(z, 10£); or 

(ii) B is extended (D, K) -regular. 

Let V be any Dirac operator with magnetic field B. Set P := e~ 5 £~ 2 , R[P) '■— (V 2 + P)~ x , let 
< ip < 1 be a function with dist(zo, supp if) > 2£. If e < s{K), then the following estimates 
hold for any u G D 

tr R[P] 2 (u,u) < c(\B(u)\P- 3 / 2 + P- 1 ' 2 ) (7.5) 
tr(R[P}V^p 2 VR[P]\{u,u) < c(\B(u)\p- 1/2 + P 1/2 ) , (7.6) 

where recall that tr := trca stands for the trace in the spin space. 

Remark. The diagonal elements in (|7.5|) . (|7.6|) are gauge invariant, i.e, they do not depend 
on the choice of the vector potential in the Dirac operator. 

Using these Propositions, we can complete the proof of the estimate (|4.11|) in Theorem 
14.31 We use that the function F(x) := P(x) and the operators A := V, Ai := T>i satisfy the 
conditions of Proposition 17.11 by using ()5.9|) . Setting n = E in ()7.4|) we obtain 

u^V^Rp^Rp+eV^^E^&E (7.7) 

2 n (cV 1 ' 2 g P 2 9 2 (pr'R^V^V^] + R 2 [P^ dJV 1 ' 2 , E^j dE 

^Tr VeffPiRiiPilVupfDiRiiPil + P?Rl[Pi] 



< 



c 

i<El 



using that J °° n(T, E)dE = Tr T for any positive operator T. This sum of traces can be 
estimated by 

cJ2 I V(x)9t(x)(\B l (x)\Pl /2 + Pf^dx < c I VP l / 2 {\B\ + P) (7.8) 

using Proposition 17.21 with D = Di, £ = £i, B = Bj, K = 100 and for sufficiently small e. The 
construction of B« for both weak and strong indices in Section IS31 guarantees that either (i) or 
(ii) holds true in Proposition 17.21 We also used that |Bj| < c(|B| + Pi) and Pi < cP(x) on the 
support of 9i (see (15. 1)1 . ()5.3|) and (j5.7|0 . moreover that Y^i®t — = 1- This completes 

the proof of (ETTTJ) . 
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8 Cylindrical geometry of the second localization 



Throughout this section B = (Bi, B 2 , B3) is an extended (D, X)-regular magnetic field. Let 
Bqc be the value of B outside D, we set b := |Boo| and := Boo/6. The corresponding 
magnetic 2-form (3 is given by 

(3 := B 3 dxi A dx 2 + B 1 dx 2 A dx 3 + B 2 dx 3 A d^i . (8.1) 

Let Zq be the center of D and define the supporting plane of B, 

V := {z G R 3 : (z - Zq) ■ noo = -£} , 

to be the plane that is orthogonal to the parallel field lines outside D. We fix an ortho normal 
basis pi,P2 in V such that pi,p 2 , *ioo is positively oriented. Any point z in V can be identified 
with a point z = (z\, z 2 ) G R 2 via z — z Q = z\p\ + z 2 p 2 , i.e. Zi = pi ■ (z — Zq), i = 1, 2. We will 
use these coordinates to parametrize V . 

For any z G V we denote by v 9 z( r ) the field li ne through z with arc length parametrization 
r, i.e. 

<fz(r) = -^fzir) = n[<p z (r)) , <p z {0) = z . 

Since B is extended (D, i^)-regular, ip z is constant outside D. If e is small enough (depending 
on K), we can assume that the length of the field line within D is at most 4£ using ()5.2j) . 
Therefore 0z( r ) is constant for \t\ > Al. 

Every field line intersects V since n nowhere vanishes and n • > 1 — ||n — n M || 2 > | if e 
is sufficiently small. Therefore the field lines {<£> 2 (r) : z G V} form a foliation of R 3 and for 
each x G R 3 we denote by ir(x) &V the unique point such that x = yj 7r ( :r )(r) for some r G R. 

8.1 Coordinates and conformal factor 

The following lemma will be used to introduce coordinates, £ = (£1,^2, £3), on R 3 associated 
with the field line <p z { T ), z G V, which will be called the central field line. The field line will 
be characterized by £1 = £2 = 0. The point z G V will be called the base of the coordinate 
system. The coordinates are functions of x G R 3 and the inverse function will be denoted by 
x(£) : R 3 — > R 3 . We may also use the notation £ z (x) and x z (£) to indicate the dependence 
on the base. For notational convenience we sometimes use £j_ := (£i,£ 2 ). 

In order to treat different error terms we introduce a notation similar to the standard 
"big-oh" notation. 
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Definition 8.1. Let k and a be nonnegative integers and let £ > be a real number. We say 
that a complex function /(£) is of class C^(|£j_| a ) if there exists a constant C such that 

|<9 e ra /(OI <crH 

for any multiindex m = (mi,m2,m^) with |m| := mi + ni2 + < k. The definition can 
be extended to matrix valued functions and to forms by replacing the absolute value with any 
matrix or form norm on the left hand side of llS.fy) . For £ = 1 we set := C^(|£±| a ), 

fork = we set O e (\£ ± \ a ) := O e (\£ ± \ a ) and for a = we set 0{{l) := O£(|£ ± | ). 

Remark. With a slight abuse of notation Ofcd^jJ ) will be used to denote not only the 
class of these functions but any element of this class, similarly to the way 0(1) is used. 

We also would like the magnetic field to be of constant strength along the central field line 
which is achieved by a conformal change of metric with a factor Q. Let ds 2 be the standard 
Euclidean metric and ds^ := Q 2 ds 2 be a conformally equivalent one. The following lemma 
describes the necessary information about the new metric and coordinates. The proof is given 
in Section liOl 

Lemma 8.2 (New metric and coordinates). Given positive numbers K,£, a ball D of 
radius £, and center z , an extended (D, K)-regular magnetic field B ; an orthonormal basis 
Pi,P2 in the supporting plane V such that pi,p 2 ,n 00 is positively oriented and the coordinate 
identification z G V <-> z = {p\ ■ (z — z ),p 2 ■ (z — z )) G R 2 . If e small enough depending on 
K, then for any zfP there exist coordinate functions £ = £ 2 (x) = (£1,^2, £3) = (£1., £3), an d 
positive functions h(£) G C 2 (R 3 ) with the following properties: 

£*(x) = and £(z) = x- z for x G V , (8.3) 

ei(^(r)) = Vr . (8.4) 

The function (z,x) G R 5 1— > £, z {x) G R 3 and the inverse function (5, £) — > x 2 (£) belong to 
C 3 (R 5 ). Moreover, if and Dx denote the Jacobians of these functions, then for 7 = 1,2,3 
we have 

\\D^x\\, ||L> 7 £|| < c(K)e 7 r 7+1 . (8.5) 
The metric ds^ := Q 2 ds 2 can be expressed as 

2 

d4 = "o^Kj + ^ 2 def , (8.6) 



mm 



(8.2) 
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where e C 2 (R 3 ) satisfies 

sup {|oy(0 - <%| : £ < |al < 10^} < c(K)e . (8.7) 
Moreover, = 5^ away from the set £ < < 10£ ; i.e. 

ds 2 = d£ 2 + d£ 2 + /i 2 d£ 2 on the domain \£ x \ < £ or |£J > 101 . (8.8) 
The functions Q, and h also satisfy 

Q = a.ij = h = 1 on the domain \^\ > 3£ , (8-9) 
n = m)(l + eO«(\!; ± \)), and h = 1 + eOi{\^\) 



with 
and 



|B(x(0,£ 3 ))| 
b 



1/2 



8.10) 
8.11) 
8.12) 



£1 = 1, d^ ± h = 0, on the domain > 10£ . 

Globally, the following bounds hold 

||/i-l|U, ||n -1 ||oo <c(iT)e, (8.13) 

IIVX-IU, IIV^IU, ||V 7 fi||oo < c(^r , 7 = 1,2. (8.14) 

Moreover, there exists an orthonormal basis {ei, e2, 63} in i/ie ds^ metric such that e% = h~ l d^ 3 
everywhere, and ej = d^, j = 1, 2, apart from the region {£ : £ < < 10£, |£ 3 | < 4£}. 

Remark. The estimates in Lemma 18.21 actually depend only on the variational lengthscale 
of the direction of the magnetic field, n, and they are independent of the variational lengthscale 
of its strength. The proof given in Section lA. 41 uses a construction that relies only on the field 
line structure and on the logarithmic gradient of |B| along the field line. However, 

— div n 



V 


n | 


B| 




B 





since B = |B|n is divergence-free, therefore derivatives of n alone control the errors. 
We define the spin-up projection associated with a field line through a given point. 
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Definition 8.3. Given z G V, the field line tp z (r), the associated coordinates £, z (x) and the 
inverse function as defined in Lemma \8.&A Then the spin-up projection associated 

with <Pz(t) is given by a 2 by 2 matrix 



Pl(x):=±[l + <T-n(x z (0,tl(x)))] (8.15) 



at any point x G R 3 . 



Note that P] is constant on the level sets of £ 



8.2 Cylindrical partition of unity and grid of field lines 

We start with a technical lemma. 

Lemma 8.4. Given y G V and the associated coordinates {£,1}, k = 1,2,3, as constructed in 
Lemma [8~2] then for any sufficiently small e < s(K) and any z G R 3 

1 - <^i%<2, (8,6) 

2- \y-n(z)\- 1 ; 

\\P}(z) - Pl {z) {z)\\ < cKer'ly - ir(z)\ , (8.17) 
where \\ ■ \\ denotes the standard norm of 2 by 2 matrices. 

Proof. Denote u = ir(z) and set q{r) := £ y ((p u {r)) — £ y ((p y {j)) e R 3 an d ^ r ij) '■= 1±( r ) — 
(9i( r )> 0.2(7)). We have |g(0)| = |r(0)| = \u — y\ and by <f(r) = n(<^(r)) and Lemma 1^21 we 
can estimate 

\q(r)\ < \\D x aoo\Vu(r)-Vy(r)\ + \\D 2 MUMr)-Vy(r)\ 

< (HA.eiUlVnlU + \\D 2 x C\\oo)\Mr) - 9y{r)\ 

< (ii^eiuiivniu + n^eiioo)ii(^e) -1 iiock(r)i 

< c(K)er x \q{r)\ (8.18) 

for |r| < 4£ and q(r) = for |r| > 4£. Therefore sup T |?(r)| < \u — y\e c ^ K ' £ by Gromwall's 
inequality and 

sup |r(r)| < sup \q(r)\ < cKet~ x e c{K)e \u - y\ . 

T T 

Combining this with |r(0)| = \u — y\ we obtain \\u — y\ < sup T \r(r)\ < 2\u — y\ if e is 
sufficiently small. Note that for some r 

|ft(z)|= el(^)(r))-ei(^(r)) =|r(r)|, 
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which concludes the proof of (|8.16|) . 

For the proof of (|8.17|) we again set u = ir(z) and by Definition 18.31 and Lemma 18.21 we 
estimate 



\\puz)-pi(z) 



< 



iVnUoo 
|Vn|L 



< iivniu(|^(o,ef(^))-^(o,e 3 ») + xy(o,e 3 (z))-x u (o,Q(z)) 

< ||Vn|| 0O (2||^riU + ||A^1U)|y - u\ 

< cKe£^ 1 \y — u\ , 

using \x y (t:) — x y (^')\ < 2|£ — £'| that follows from (J8.5)) if e is sufficiently small. This completes 
the proof of (IHT7|l . □ 

We construct a grid of field lines. Choose a square lattice y := : j 6 Z 2 } on P 
with spacing 6 _1//2 , i.e., \yj — yk\ = b~ l l 2 \j — k\, j, k E 7? . Applying Lemma l%^2l to each field 
line tpy., we construct conformal factors fij, orthonormal bases {e^\ e% } and coordinate 
functions = $\ t?)- We 

now construct a set of Gaussian localization functions 
with a transversal lengthscale of order b" 1 ^ 2 that are essentially supported around the field 
lines ifyy Let rj < \ be a small positive number to be specified later and we define 



Vj (x) = exp ( - ^ { l\x)f 



(8.19) 



We set Pj(x) := Py-(x) to be the 2 by 2 spin-up projection matrix associated with the field 
line through yj (see Definition I8.3j) . 

Lemma 8.5. If e is sufficiently small depending only on K, then for any 7 > ; K > we 
have 



(8.20) 
(8.21) 



uniformly in x G R 3 . Moreover, there is a universal constant Cq and for any < A < 1 there 
exists < 77(A) < I such that for any rj < 77(A) 



iGZ2 



b[X-v Z b[^i\x)Y)PUx)+Q 



0-2 



> 



(8.22) 
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Proof. Since (?7&£ 2 _) K exp(— < 0(7, k) uniformly in it is sufficient to estimate vj 
for the proof of (|8.20j) . Using ()8.16|) we obtain 

J2 v] /2 (x) < J^exp ( - - vr(x)| 2 ) < c(rry) -1 

since runs through a square grid with spacing 6 -1 / 2 . The proof of (|8.21j) is similar. 

For the proof of ()8.22)1 we define k G Z 2 to be an index such that \yt~ 7r(a;) | < fe -1 ^ 2 . Then 
by ()8.16|) and Schwarz' inequality 

|y fc - %f < 2&- 1 + 2|vr(x) - y,f < 26" 1 + 4|£i(x)| 2 . 

Combining this estimate with (|8.17|) and using that (P^f = P^ we have 

p] > \pI - 2\\p] - p\f > \pI - r*\y k - Vj \* > \pl - 4 r 2 |£il 2 - cb-n- 2 

and 

P) < 2Pl + 2(Pj - Plf < 2Pl + 4r 2 |£i| 2 + c&- x r 2 

if e is sufficiently small. We omitted the x argument for brevity. Therefore we can use (|8.20j) 
and (|8.21|) to estimate 

> \ E<( A - v&i£iM - cj2^(\b\e ± \ 2 + c \+ v ^{\ 4 + ^ieii 2 )r 2 

> h _^X-c)Pl-c V - l t- 2 

if 1] is sufficiently small. We can choose Cq to be the universal constant c in the last formula 
and the proof of ()8.22j) is completed. □ 



9 Dirac operator on R 3 with a general metric 

The following sections summarize basic information about the Dirac operator over a non-flat 
manifold. More details are found in [ES-IIIj (the sign of A is chosen to be the opposite in 
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this paper). The presentation here is simplified because the spinor bundle is trivial and we 
can work in a global orthonormal basis. 

Throughout this section we shall consider R 3 with a general Riemannian metric g = (-, •) 
and we shall consider the Dirac operator for this particular Riemannian manifold. The Dirac 
operator will be an unbounded self-adjoint operator in L 2 (R 3 ) ® C 2 (the subscript g refers to 
the fact that the measure is the volume form of g). 

Let {ei, e 2 , e 3 } be a global orthonormal basis of vectorfields and let {e 1 , e 2 , e 3 } be the dual 
basis. If X is a vectorfield on R 3 , we denote by 

1 r 3 i 

3=1 

the spin projection in direction X with respect to the basis {ei, e 2 , e 3 }. 
We also introduce a covariant derivative on L 2 (R 3 ) ® C 2 by 

Vx := d x + | <r • . (9.2) 

Here we define 

w(X) := ((Vxe 3 ,e 2 ), (V x ei,e 3 ), (V x e 2 ,ei)) , (9.3) 

where Vx refers to the Levi-Civita connection on vectorfields for the metric g on R 3 . 

If a is a (real) 1-form we define the corresponding covariant derivative on L 2 (R 3 ) ® C 2 
(see Proposition 2.9 in [ES-III ) 

V" :=V x + ia(X). (9.4) 
The magnetic 2-form is (3 := da. We define the Dirac operator by 

3 

P a := J>'(-iVy ■ (9.5) 

3=1 

It is a symmetric operator in L 2 (R 3 ) ® C 2 (Theorem 3.2 in ES-III]). Note that P a also 
depends on the metric g and the choice of {ei, e 2 , e 3 } but this fact will usually be suppressed 
in the notation. 

For notational convenience we introduce the following vector of covariant derivatives 

IT :=HV£,-^ 2 ,-zV£). (9.6) 
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With this notation we may write V a = a ■ IT*. Note that the components of 11° are not 
self- adjoint, however the components of the vector 

D a = (D« D£, D«) : = IT - % - (div ei, div e 2 , div e 3 ) (9.7) 

are self-adjoint operators. 

For any one form A = Aie 1 + A 2 e 2 + Ase 3 we define cr(A) := Aicr 1 + A 2 cr 2 + A3O" 3 . The 
Lichnerowicz' formula (see, e.g., Theorem 3.4 in |ES-IIIj ) states that 

[v a f = [ir]* -n a + -r + <t{*/3) , (9.8) 

where R is the scalar curvature of g and * denotes the Hodge dual. 
In terms of operators, the Lichnerowicz' formula reads as 

3 3 

[V a f = [Df + \R+\ £>v e / + ^ ( div e i) + °W) • (9-9) 

3=1 j=i 



For the flat Euclidean metric with the standard ortho normal basis we have u> = 0. In this 
case if a denotes the 1-form dual to the vector field A = (Ai, A 2 , A 3 ) then II a = p A where 
is the vector of operators {—id\ + Ai, —id 2 + A 2 , —id 3 + A 3 ). Therefore we obtain the usual 
Dirac operator T> = D a defined in Section Q with (3 = da given by flEU). Moreover, P x = P] 
if X(x) = n(x 2 (0, £ z (x))) from (pTTKj) and JSHJ). 



9.1 Gauge transformation 

In the previous construction V a and n a depend on a and also on {ex, e 2 , e 3 }. Up to a unitary 
equivalent gauge transformation, however, T> a and n a depend only on the metric g and the 
magnetic 2-form (3. Similarly, the spin projection Px defined in (19. lj) is gauge- invariant. 

More precisely, given another 1-form a 1 with da' = (3 and another orthonormal basis 
{e[, e' 2 , eg} with the same orientation, we denote the corresponding operators by T>' and IT' and 
let P' x be the spin projection. There exist a real valued function <p{x) and a continuous function 
R(x) G 5*0(3) on R 3 such that a' = a + dcj) and J2k w k e 'k = X]fc(-^ w )fc e /c f° r an Y w ^ R 3 - 
Let Ur{x) G SU{2) denote the image of R(x) under the isomorphism 5*0(3) — > 577(2)/{±l}. 
The requirement that Ur(x) be a continuous function of x determines Ur uniquely up to a 
global sign. In particular 

U R {a-v)U R = cr-{R-v) (9.10) 



34 



for any v G R 3 , i.e. R(if),eTif>) = (Uri/;, ctUriP) for any if) G C 2 , where (if>,crif)) denotes the 
vector (^(ipja 1 ^), (ip,a 2 ip), (ip,a 3 ip)^ G R 3 . 
We define the unitary operator of the form 

[Ur^x) = e^ x) U R (x)iP(x) , (9.11) 

then 

P'x=Ur^PxU r ^ (9.12) 

and 

V' = U R4 V a U R ^, and w II' = (Rw) ■ Wr^TFUrj (9.13) 
for any w G R 3 . In particular, the spectrum of T> a and the functions 

tr f -r— — — (2,0;)^ and trf — — V a ip 2 V a ——^- (x.x 

depend only on g and f3, where <p is any function on R 3 and c > is a constant. 

9.2 Change of the Dirac operator under a conformal change of the 
metric 

Let Obea positive real function on R 3 and let gn := Q 2 g be a metric which is conformal to g. 
Consider the (fi, f2, fs) '■= (^ _1 ei, Q~ 1 e2, Q^e^) orthonormal basis in g^. Given a 1-form a 
we let and Vq denote the corresponding covariant derivative and Dirac operator. With 
the notation 

EE := ( " *VJ°, -iVg n , -iV% n ) (9.14) 
we have = cr ■ II£. Then from Section 4 of |ES-lllj 

V£ = tt- 2 V a tt (9.15) 

and 

W a f = V a x + iCr>pf*), a(dQ)} (9.16) 

for any vector X, where X* refers to the 1-form which is dual to the vector X relative to the 
metric g, and a(X*), a(dfl) are computed in the {ei, 62,63} basis. In particular, 

II£ = Q- 1 ^ - {Q- 2 ([a 1 , a(<m)}, [a 2 , <r(dfi)], [<x 3 , <r(dfi)]) . (9.17) 
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9.3 Constant approximation of the magnetic field along a field line 



The goal of this section is to express the Dirac operator with a non-homogeneous regular mag- 
netic field as a sum of a constant field Dirac operator and some error terms in a neighborhood 
of a given field line. This can be done if the original Dirac operator is already written in an 
appropriate orthonormal basis and with a carefully selected vector potential. The basis and 
the vector potential are determined by the local magnetic field. 

Given an extended (D, i"T)-regular field B. Consider the corresponding 2-form f3 and a 
fixed field line. Let the coordinates (£1,^2, £3), the new metric = fi 2 ds 2 with a conformal 
factor Q and the orthonormal basis {ei,e2,es} be as constructed in Lemma I5~2l associated 
with the given field line. Let a denote a vector potential, da = /3, to be chosen later. Let n a 
be given by (JUSD and let V a := cr • IT*. 

On the central line and in the regime > 10£ the magnetic field f3 is constant in the 
dsft metric: 

(3{e u e 2 ) = n~ 2 (3(n ei ,ne 2 ) =n- 2 \B\ = b, P(ej,e 3 ) = 0, J = 1,2. 
This observation gives rise to the following definition. 

Definition 9.1. Given afield line, the associated coordinate system £ and the conformal factor 
Q as above such that magnetic field (3 is constant in the ds^ metric with strength b = /3(ei, e 2 )- 
Then the magnetic field (3 C given by 

/3 c :=&d£ 1 Ad& 

is called the approximating constant magnetic field along the field line. 

The magnetic field /3 C is clearly constant in the d£ 2 = Y^=i d£j metric. A convenient 
gauge is defined as a c := f [£id£ 2 — ^d^i], then (3 C = da c . 

In particular j3 = j3 c along the central line and in the regime > 10£. We compute the 
norm of the difference field 5/3 := (3 — j3 c and the norm of its derivative in the d£ 2 metric. 
Using JHHJ), dHZZJ), (EH), (EH and (fPJl . (jfTKjl we obtain 

6(3(0 = ebOi(\^\) (9-18) 

and 5/3(0 = if > 101 

Next, we define an appropriate gauge a for the original magnetic field /?, da = /3, such 
that a — a c be small. The following Lemma was given in jES-1] (Proposition 2.3). Although 
it was stated in a slightly weaker form, the explicit formula (2.20) of |EH-Ij gives the following 
stronger result with a straightforward computation: 
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Lemma 9.2 (A-formula). Given any C 2 magnetic 2-form (3 on R 3 with Euclidean cooordi- 
nates (^1,^2,^3)- For k,m e N we define 



bk,m(£,l 



u k sup \\V m [3(u,z 2 ,z 3 )\\du + / u fc sup ||V m /3(2 1 ,M,2 3 )||dM 

II 2:2,^3 JO *1,«3 



JTien t/iere exists a 1-form a generating (3, da = (3, such that 
||a(OII < c[6o,o(a)+6i,i(^ , 
||Va(OH < cfsup{||/3(u)|| : < + M£i) + 



(9.19) 



(9.20) 



We apply this lemma to the magnetic 2-form 5(3 and we denote by 5a the generating 
1-form. We define a := a c + 5a, then a generates the original magnetic field (3, da = (3 and 
it is close to the linear gauge a c of the constant field (3 C using (J9.18J) and Lemma I9~2l 



(a-a c )(0 = eb£O f 1 (\^\ 2 ). (9.21) 
The norm of the left hand side is computed with respect to the standard metric. 

Definition 9.3. With the notations above, the Dirac operator 

3 

V:=J2A~^ k +a c (da (9.22) 
fc=i 

with a constant field (3 C in the d£ 2 metric will be called the approximating constant field Dirac 
operator along the field line. 

By the properties of the coordinate vectorfields d^ k and the orthonormal basis {ei, e 2 , e 3 } in 
the gn metric from Lemma IHT21 and by the definitions (|9.2|) . (|9.4|) . ()9.6|) we have, for sufficiently 
small e, 

3 3 
V a = a- ■ IF = V + A^k) - a c (9&)] + ^H%) + M o , (9- 23 ) 



fc=i 



where /C = (/Ci, /C2, /C3) and K.f-, M.o are 2 by 2 matrix valued functions. We use the bounds 
(IQl . (EHUL dHIIH and the estimate (I9~2T1) to obtain, for k = 1, 2, 3, 



[a(e fc ) -a c (%)](0 



(9.24) 
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We obtain from flQSj) and flOU) that 

3 

V a = V + ^/C,H%) +.M (9.25) 

k=l 

with matrix valued functions that satisfy 

M = b£0{(\^\)+ r^(l), /C 1)2 = <3(|£lP) , ^3 = ^(1^1) (9.26) 

for any 7 > if e < e(K). These estimates follow from Lemma [8 .2[ especially from the fact 
that efc = <%, k = 1,2, apart from the region £ < < 10£, where ()8.6|) holds, i.e. /Ci,/C2 
are supported in this region. 



10 Positive energy regime: Proof of Proposition 16.2 



We first notice that both sides of ()6.5|) scale as £~ 2 , hence it is sufficient to prove the result for 
i := 1. We can apply the constructions of Section |H] for the magnetic field B to obtain confor- 
mal factors Qj, orthonormal bases {e^, e% coordinate functions £^ = ,£2 ,63 ); 
spin-up projections Pj and Gaussian localization functions Vj concentrated along the field line 
passing through yj. We recall that yj was a lattice with spacing b~ x l 2 on the supporting plane 
(see Section IS^j) . 
We first estimate 



£ -^[ P i + CT • b] > id 2 + 



since b = \Boo\ > e 2 an d SU P |B| < 2b if e is sufficiently small. 



Using this estimate and (j8.20J) - (j8.21|) we have 

V 2 + us- 5 - M X 2 V > l(v 2 + e- 2 lp 2 A + fie' 5 -2M X 2 V 

> cr]J2 ( Vv j V + £ ~ 2 |pa ■ v 2 p A + fie~ 5 Vj - cM X 2 Vv\ 



(10.1) 



if e is sufficiently small (depending on /1). Notice from the explicit formula (J8.19|) that 





2 


\PA,Vj] 


= \V Vj \ 



(10.2) 
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Therefore by Schwarz' inequality pa • v ]pa > \ v jPA v j ~~ ci]bvj, and using this estimate in 
IjlO.lJ) . including the negative error term into ne~ 5 Vj and subtracting the pointwise inequality 
(|8.22|) we obtain for any < A < 1, rj < 77(A) (see Lemma [8.5(1 that 

V 2 + fie- 5 - M X 2 V (10.3) 

> cq (H V + £ ~ 2 ± b v jP 2 A v J ~ V P( X - r, 2 b[&f)P} + - cM X 2 Vvf) . 

j 

The error term CqvI- in (|8.22(1 has been absorbed into fie~ 5 Vj if e is small enough depending 
on 77, /i. The inequality (|8.22(1 has been subtracted to prepare for a later step. Hence 



Tr(£> 2 + ij£- 5 - MxV) j (10.4) 
< cr] j Tr + e-^vjplvj - vp(X - v 2 b[^l ] } 2 )P] + f £~S - cM X V^ 



by |TrQ^. < 52- |Tr(^)_| that follows from the variational principle. The following 

lemma is the cylindrically localized version of Proposition 16.21 



Proposition 10.1. With the notations above and setting W := cM\ 2 V we have 

Tr(vvp + £- 2 ^p A ^ - vjb(2- 7 - r} 2 b[ff] 2 )P] + f<r 5 ^ - Wv^ J (10.5) 

< c J ([v 2 Wf' 2 + b{v 2 Wf' 2 \ 

for each j if e is small enough depending on K,^i,r]. 

Choosing A := 2~ 7 and rj := r](2~ 7 ) (see Lemma f8. 5(1 . Proposition 16. 21 directly follows from 
this proposition, from (j!0.4() and ((8.20(1 . □ 

Proof of Proposition UUTR The proof contains three transition steps that are performed locally 
around each field line from the grid constructed in Section 18.21 First we replace T> with a 
Dirac operator £>q in a metric that is conformal to the Euclidean one. The conformal factor 
Q is chosen such that the strength of the magnetic field becomes constant along a field line. 
Second we replace the volume form dx with the volume form d£ = d£i A d^2 A d^3, where 
£ is the coordinate system associated with the chosen field line. Then we perform a gauge 
transformation so that becomes close to the Dirac operator T> := cr ■ (—id^ + a c (d^)) with 
a constant magnetic field (3 C = da c = b d£i A £ 2 in the linear gauge a c = |[£id£2 — 62^1]- 
Finally, the operator T> can be analyzed explicitly. 
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For each fixed j we consider the constructions in Section El with the metric gj '■= gn = 
Q 2 dx 2 . We shall apply Section T9.2I to the Euclidean metric with the standard basis vectors 
{di,d 2 ,d 3 }. The vectors {flj 1 di,flj 1 d2,ftj 1 d 3 } form an orthonormal basis in gj. Let T> a - and 
ITj denote the corresponding Dirac operator and the vector of derivative operators as defined 
in (J9.15|) and ()9.14|) . We recall that a is the 1-form dual to the vector potential A in the 
standard metric. Using the estimates in Lemma f8. 21 to the formula ()9.17|) . we obtain 

n; = nj l p A + e0 2 (l) , and V] = nfV + e0 2 (l) , (10.6) 

where the error terms are functions. By Schwarz' inequality we obtain the following pointwise 
bound 

|pA^| 2 >||n^r-c(ir)e 2 ^| 2 and \V^\ 2 > §|P>| 2 - c{K)e 2 \i,\ 2 (10.7) 
since | < Qj < 2 if e is sufficiently small. Therefore 

Vjj&vj > Iv^}* ■ njvj - c(K)e 2 v 2 (10.8) 

and 

Vvp > | [V]]*vp] - c{K)e 2 v] , (10.9) 

where star denotes the adjoint in the standard L 2 -space. By applying the inequalities ()10.8j) 
and (|10.9|) we have 

Trhvfi + e - 2 ^ b v jV 2 A v 3 - vp(2- 7 - r) 2 b[C±] 2 )P] + f e~\ - Wvfj | (10.10) 

< iK^iN^ + e ~ 2 i v ^T ■ n a jVj - v*b{± - Wh[&f)P] + - %Wv] 

The error terms in ()1U.8|) and ()10.9[) have been absorbed into the ^e~ 5 Vj term using Vj < 1 if 
e is sufficiently small. 

The right hand side of (jlO.lOj) is invariant under an SU(2) x U(l) gauge transformation 
Ur,<$> as defined in ()9.11|) . We shall choose R to be the rotation from {ft,J 1 di,Qj 1 d2,Qj 1 d3} 
to the basis {ei,e2,es} constructed in Lemma ET21 and to be such that a = a + d0, where 
a is constructed in Section l9~31 In particular becomes := ~[1 + cr 3 ] according to (J9.12|) 
since n = e 3 along the central field line. Therefore the right hand side of (110. 10J) continues as 

(flZHZi) = | Tr (\D a ]*v A V a + e~ 2 j- b v[n a ]* ■ U a v - v%{± - 8rj 2 b£ 2 ± )^ + ^e~ 5 v - 8Wv 6 ^j , 

(10.11) 

11 11 • QTlH h i 

3 ' 3 ' 



where we also omitted the j index for brevity, i.e. T> a = V®, U a = II", v = Vj and £j_ = ^j 7 



for the rest of this section. 
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Now we translate our problem from the standard L 2 (dx) space to the L| := L 2 (d£) space. 
We change the measure from the volume form dx to d£ = d£i A d£ 2 A d£ 3 . Since these two 
volume forms are comparable by a factor of at most 4 by (j8.13|) of Lemma l£T2l if e is sufficiently 
small, we can use Lemma fA.81 from the Appendix to obtain 



(IHHTJ) < | Tr L | ([V a ]*v 4 V a + e~ 2 ±v[Tl a }* ■ ITv - v%\ - 8r/ 2 ^i)^ T + - 32HV 

where the trace and the adjoints are computed in the L 2 (d£) space. 

We remark that already on the right hand side of (jlU.lOj) it could have been natural to 
transform the trace on L 2 (dx) to the trace on L 2 (f2 3 dx), according to Lemma fA.81 

We introduce the function 

G = G(0 :=1 + \/&min{|£ ± |,l} , (10.12) 

and we notice that 

sup G(0M0" < c(p, q)v~ p/2 , p > 0, q > . (10.13) 

We recall the definition of Df, and the decomposition (|9.25|) from Section IQ1 Since £ — 1, 
the estimates ()9.26|) are translated into 

M = 0(G 2 ), \\VM\\ = VbO(G), /C 1 , 2 = 1 (|ar), £ 3 = 0i(|£i|). (10.14) 

for any 7 > 0. Here C(G fe ) denotes the class of functions F(£) on R 3 with sup^ |F(£) |/G fe (£) < 
00. 

Hence 

J v A \V a m > \ f \v 2 V^\ 2 di - C E / \v 2 ^ k d^M ~ J v A 0(G A M 2 dt . (10.15) 



In the second term on the right hand side we first commute v through the derivative. Notice 
that [<%,t>] = r}bvO(\£±\) for k — 1,2 and [<%,f] = 0. Then we use the estimates (|10.14j) and 
(110.131) to obtain 

33 2 
J2 [ \v 2 JC k d^\ 2 d£ < 2^ / v 2 \\}C k \\ 2 \d^\ 2 d^ + 2j2 [ ^M 2 [??&0(|6_|)] 2 M 2 d£ 

k=l ^ k=l k=l J 

3 

< cb- x J2 [ l<%Wfd£ + c / v 2 \i>\ 2 d£ (10.16) 
k=i •* 
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if e is sufficiently small. From the last term in (|1(J.15|) we obtain a similar error term as in 
flflHBD using IjlU.lSj) . hence 

j v^V^&i > \ J \v 2 V^\ 2 d£ - cb- 1 J - c j v 2 \^M . (10.17) 

We also define for k — 1, 2, 3 

n„, fc := + (1 + 2??)a c (<%) 

and 

3 

V v := cr ■ Tl v = cr^n^ 
fc=i 

which is a Dirac operator with constant field (1 + 2rj)b d£i A d^2 in the d£ 2 metric. Notice that 

V v v 2 = v 2 V + 2irjbv 2 (cr% + a 2 £ 2 )a T . (10.18) 

This identity, called the magnetic localization formula, was introduced in [ES-IIj . 
Hence, using (|10.18|) to continue (|10.17|) . we obtain 

J v 4 \V a i)\&£ > \J [Dy^Ai-^rfb 2 J v 4 il\^i)\ 2 di (10.19) 

3 

- crl / (Ei%^i 2 ) d ^- c / v2 \^ 

k=l 



> 8 



3 

fe=i J J 



In the last step we used that V 2 > 26(1 + 2rj)a^ > 2ba\ i.e., that on the spin-up subspace 
{ip : a^ip = ip} the constant field Pauli operator is bounded from below by twice of the 
constant field. 

We shall control the second negative error term on the right hand side of ()10.19j) by the 
term e~ 2 (2b)~ 1 v[II a ]* ■ II a v. Notice that the following inequality holds pointwise 

|n>| 2 > i(X>^| 2 ) - 4|H0II 2 H 2 - 4(su P |M) 2 H 2 

k=l 
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using flH]|, ((231) and <pO>jl . We can estimate ||a(£)|| < c& 1/2 G(£) from (fg^jj) and the explicit 
choice of a c . We also estimate ||u?|| < cKe by Lemma 18.21 and the same lemma is used to 
estimate the transition from Ylj l^y^l 2 to Y2j I'^V'I 2 - Therefore 



|n«vf >~(£l%Vf)-c&G 



2L/,|2 



fc=l 

if e is sufficiently small, hence 



3 

^ J |ir^| 2 d£> JL | |%,^| 2 )d£- V 2 | W 2 |^| 2 de (10.20) 

fc=i 

using (gnu . 

Combining (110. 19|) and (ll().2()jl we obtain 

[V a ]*v A V a + £T 2 ^;[ir]* • U a v - v%\ - 8r] 2 b^ 2 ± )^ + ue' 5 v 2 - 32Wv 6 

> WV 2 n v 2 - 32Wv & 

o 'I 

if e is sufficiently small depending on is,rj and K. Since Tr(X*HX)- < \\X*X\\TrH-, 

Tr(^v 2 V 2 v 2 -32Wv 6 ^j \ < § Tr (v 2 - 256Wv*\ ; (10.21) 

< cj ^(l + 2r/)(256iyi; 2 ) 3/2 + (256^ 2 ) 5/2 )d^, 

where in the last step we used the Lieb-Thirring inequality for a constant magnetic field 
|LSY-IIj . This completes the proof of Proposition llO.il □ 

Remark. The reader may have found it confusing that along the proof of the positive energy 
regime we used the Birman-Schwinger principle (|4.fij) back and forth several times. It occured 
first in ()4.f)jl . ()4.9j) . then in (|fi.7j) . (jfj.lOjl . and finally implicitly in (|l().21jl . when we referred 
to the Lieb-Thirring inequality with a constant magnetic field whose proof also relies on the 
Birman-Schwinger principle. The frequent changes back to an expression on the sum of the 
negative eigenvalues were purely for the purpose of compact presentation of the intermediate 
results. It would have been possible to use only ()4.6j) and stay with the resolvent language all 
the time since all estimates done for the operators are equally valid for the resolvents. In this 
case, of course, we could not have referred directly to the result of [LSY-II on the constant 
field case, rather to the details of that proof. 
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11 Zero mode regime: Proof of Proposition 17.21 



First notice that the inequalities in Proposition 17.21 are scale invariant in powers of £; both 
sides of ()7.5|) scale like £ and both sides of f)7.6j) scale like i~ x . Therefore we can set i = 1 for 
the proof. The arguments for weak magnetic fields and for extended (D, -fT)-regular fields are 
different. 



11.1 Weak magnetic field 

For weak fields (|7.5|) and ()7.6|) will be estimated by a universal constant c if e < s(K). We 
need the following lemma: 

Lemma 11.1. Let X, Y be self-adjoint operators such that X > 0, X + Y > and \\Y\\ < M 
for some constant M > 0. Then 



1 



< 



< 



(X + Y + 2M) 2 ~ (X + M) 2 - X 2 + M 2 
Proof. Consider the resolvent expansion 

1 1 1 



11.1) 



X + Y + 2M X + M X + M 
hence by Schwarz' inequality 

1 2 1 



(Y + M) 



X + Y + 2M 



(X + Y + 2M)< 



< 



< 



< 



+ 2- 



(X + M) 2 X + M 
2(2M\ 



(Y + M) 



(X + Y + 2M) 



(X + M) 2 
4 



X + M 



(Y + M)(Y + M) 



:(Y + M) 
1 

X + M 



1 



X + M 



(X + M) 2 

since (Y + M) 2 < 2\\Y\\ 2 + 2M 2 < AM 2 . By positivity of X we have (X + M) 2 > X 2 + M 2 
which completes the proof. □ 

If || B || oo < ce~ 2 then using Lemma fl 1.11 we obtain 



R[P} 2 (u,u) 



'—iV + A) + er -B + e 



-5 



-2 



(u,u) < 4 



'— iV + A 



i 2 + K 5 



-2 



(u,u) 



if e is sufficiently small. By the diamagnetic inequality we can continue this estimate as 



R[P} 2 (u,u) < 4 



A + |e- 5 



-2 



(u, u) < ce . 
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This proves (|7.5|) . 

For the proof of (|7.6[) we define a smooth function < x < 1 such that x{ u ) = 1j | Vxl — c ) 
|V 2 x| < c and supp(x) H supp((/?) = 0. For brevity we set i? := -R[P] with P = e~ b . Since 
the inequality (|7.6|) is gauge invariant, we can choose the Poincare gauge A centered at z to 
generate the magnetic field. In particular, IjAHoo < ce~ 2 since by assumption B is supported 
on D and ||B||oo < ce~ 2 . 

Let {X, Y} := XY + YX denote the anticommutator. Notice that [R, x) = R{T>, [V, x]}R 
and 

{V, [V, x)} = (-i){<r ■ HV), a ■ V X } + (~i){<T • A,cr ■ V%)} . 
We can compactly write 

3 

{P, [V, X ]} = $^(-^V + A.)jK,j + M , 
i=i 

where the JCj and are 2 by 2 matrix valued functions and from the estimate on A and the 
derivatives of x we easily obtain that sup x ||/Cj(x)||, ||A^(x)|| < ce~ 2 . 

Therefore commuting x through first, estimating if 2 < 1, then using RT> 2 R < R and 
applying a Schwarz' inequality we get 

xRVv 2 VRx = R{V,[V,x\rRV^ 2 VR{V,[V,x}}R 

< R{V,[V,x\VR{V,[V,x\}R 

3 

< ARM*RMR + 4 ^ RK.*(-iV + A)jR(-iV + A)jK.jR . (11.2) 

In the first term we use R < P^ 1 = e 5 for the middle resolvent then we use the boundedness 
of M. to arrive at the resolvent square, R 2 , that was estimated above in the proof of ()7.5|K 

In the second term we use that V 2 + P = (— iV + A) 2 + cr • B + P > (— zV + A) 2 + since 
||B|| < ce~ 2 < f = if £ is sufficiently small, therefore (— iV + A)jR(— zV + A)j < 1. The 
estimate of the second term then can be completed by using ||/C*/Cj|| < ce~ 4 and referring to 
the estimate of the square of the resolvent in (J7.5j) . This finishes the proof of Proposition 17.21 
for the case of weak magnetic field (case (i)). 

11.2 Strong magnetic field 

Here we prove Proposition 17. 21 for the case (ii). Throughout the proof we fix u and let z = n(u) 
be its base point on the supporting plane V. Consider the construction of Lemma 18.21 in 
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particular the coordinate functions £ = {£1,^2, £3) an d the conformal factor Q. We know that 
~ < Q < 2 and ||Vfi||oo < 1 if e is sufficiently small. Recall that we set t = 1, therefore the 
bounds on the right hand side of ()7.5|) and ()7.6|) saturate to c(e, K)b since P = e~ b and |B(u)| 
is comparable with b := |Boo| > 1. 



11.2.1 Transformation into good coordinates 

Similarly to the positive energy regime in Section El we perform three transition steps. We 
first change T> into 2^ := Q~ 2 UQ and the underlying measure to dx to Q 3 dx, then we change 
the measure from Q 3 dx to d£ and finally we perform a gauge transformation. Recall that T)^ 
is self-adjoint on L 2 {ds 2 ^) (see |ES-111| ). We set R^[P] := ([V^] 2 + P)- 1 and we assume that 
e is sufficiently small so that P > 2 9 . 

Then Lemma IA.9I from Appendix IA.6I states that 



tr R[P} 2 (u,u) < 2 9 tr (i$[P])£ a(n) (u,u) (11.3) 



tr \ R[P]Vip 2 VR{P])(u,u) 
< 2 12 tr (RUP}V& 2 V a n R a n [P}) L2{n) (u,u) + 2 12 Ptr (P£[P])! 2(n) («, «) (11.4) 

where the operator kernels on the right hand side are computed in the L 2 (Q) := L 2 (Q 3 dx)®C 2 
space. 

Gauge transformation of the form (|9.11|) leaves the diagonal elements of operator kernels 
invariant hence we can use the basis {ei, e2, 63} constructed in Lemma and vector potential 
a constructed in Section l9~3l to express the right hand sides of (jll.3j) . (jll.4j) using instead 
of T>q as in the proof of Proposition 110.11 

We recall that in this gauge the decomposition (j9.25J) holds, i.e., 

3 

V^ = V + J2^ k +M (11.5) 

fc=i 

with V given in (j9.22j) and K, k , M satisfy the estimates (jlU.14j) with £ = 1. 

We apply Lemma IA.7I from Appendix IA.5I to compare operator kernels on the measure 
spaces with volume forms d/z := Q 3 dx and dv := d£ = d£i A d^2 A d^. Since these two volume 
forms are comparable at every point by Lemma f8. 21 we obtain from (jA.33|) that 



ta(i$[P])Wu,u) < ctr(P^[P]*P°[P]) L2(d0 ( 



u. u I 



tr ^[P]V^ 2 V^[P]) l2 ^(u,u) < ctr ( [pV«R«{P}\ * cpV%R%[pj) ^(u, u 



L2(dO 
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where the adjoints and the operator kernels on the right hand sides are computed in the 
L 2 (d£) <8> C 2 space. 

Therefore case (ii) of Proposition 17.21 has been reduced to proving that for e < e(K) 



tr 



R&PYR&P] (u,u) < cb, (11.6) 



tr 



&VXlP]y V V«R°iP] (u,u) < cb, (11.7) 



where the adjoints and the operator kernels are computed on L 2 (d£) ® C 2 . 

Remark: The operator X>£> in general is not self-adjoint in L 2 (d£) <g) C 2 , but T> is. 

11.2.2 Proof of (HEBD and (fTT77H 

In this proof we will omit Q and a from the notation of T>^ and we will simply use T> for 
this operator. This should not be confused with the notation T> (see ()1.2|) with h = 1) used 
elsewhere in the paper. 

We recall that T> denotes the Dirac operator with a constant field (see (J9.22)) ). We also 
need the notations II = (III, n 2 , n 3 ), II = (fii, U 2 , n 3 ) from V = a ■ II, V = a ■ fl. We note 
that V can be decomposed as 

V = V ± + N (11.8) 

with 

2 

V± := ^2<T j H-j and := a 3 n 3 = (J 3 (-z%) . (11.9) 

3=1 

These operators are self-adjoint on L 2 (d£ 2 ) and 

[Pj.,n 3 ]=0, {V ± ,N} = 0, and V 2 = V>\ + iV 2 = T>\ + fl 2 . (11.10) 
We need two decompositions of the error term £ := D — T> as in (|9.25J) 

s = n-tc + M (ii-ii) 

£ = U-Z + M . (11.12) 

Here /C = (/Ci, /C2, /C3) is a vector of 2 by 2 matrices and II • JC := Xlj=i Hi^i- T ne matrices 
/Cj, satisfy p0.14|) with i = 1 and the same estimates hold for Kj, Ai as well. The following 
estimates are straightforward from (jl(J.14j) and (|9.21j) if e < e(if) 

tn is x: fc ] = o(i), [n 3 ,M} = b^ 2 o(G), j,k = 1,2,3, (11.13) 
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and the same holds for K, and M. In particular the following relations also hold: 

£ = IC-U + Mo , S = £-n + M (11.14) 

with Mo, M = 0{G 2 ). 

The following lemma collects various operator inequalities related to the diamagnetic in- 
equality. The proof is postponed to Section IX"T1 

Lemma 11.2. With the notations above we have the following operator inequalities in the 
space L 2 (dO <g> C 2 if e < e(K): 

m(— - — Vf— =- — )il < cb, (n.15) 



1 \ * 1 

1 x ! 



£ ivhr)"w^p £ £ ° [G>) (1118 > 

For t/ie constant field operator we have 

EL^- n,<c6, j = 1,2,3. (11.20) 

3 V 2 + P ~ 

The next lemma estimates the diagonal elements of explicitly computable operators with 
a constant magnetic field. The proof is given in Section fA.81 We set H(£) := min{|£j_|, 1} 
and we recall that G = 1 + VbH. 

Lemma 11.3. With the notations above and for any constants P > 1 ; b > 1 with P < cb we 
have 

tr -= (u, u) <cb . (11-21) 

(V 2 + P) 2 

For any k = 1, 2 . . m = 1, 2, . . . and /or any 2 by 2 matrix valued function T with = 
0(G(x)) we also have 

tr— !: G 2m ^ (u,u) < c(m)b, (11.22) 

V 2 + P V 2 + P 
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tr „ V [F k ]*^ T k ~ V (u, u) < c(k)b, (11.23) 

V 2 + P ^l+P 7)2 + P ~ 

tr = V H 2m ~ V (u, u) < c(m)b l - m , (11.24) 
V 2 + P V 2 + P ~ 

tr „ Uj H 2m ~ Uj (u,u) < c(m)b 2 - m , j = 1,2,3. (11.25) 
V 2 + P v 2 + P 

Let W denote either the identity I, or V, or n 3 , then 

tr = P [T k Y J^- G 2m „ W T k „ V (u, u) < c(k, m)b (11.26) 
V 2 + P V 2 + P V 2 + P V 2 + P 

(we recall that [n 3 , V] = 0). Let U denote either the identity I orV orUj,j = 1, 2, 3, and let 
< (p < 1 be a function with dist{u, supp(ip)) > 1, then 

tr = - y; 2 ^ - K^) < ce"^ , (11.27) 
£> 2 + P X? 2 + P 

tr ^ ^ lT fe r ^ - G m y 2 G m ^ J* = P < c(k,m)e- cVh . (11.28) 
£> 2 + P £> 2 + P P 2 + P P 2 + P 

Armed with these lemmas, we complete the proof of (jll.6j) and ()11.7|) . We start with 
(|11.6j) . We introduce the notation (• • -)*A for A* A if A is a long expression. All adjoints are 
computed in the L 2 space. 

We use T> 2 = T> 2 + T>£ + £T> in the following resolvent expansion: 



with 



(-)' 'z^ 3 (^ + < B > + < c >) 



V 2 + P 



ID) : ( ••• ) — X>£- ' 



X? 2 + P p2 + p 



£> 2 + p v 2 + p' 

Term (A) is explicit from (jll.21|) . In term (B) we first use f)11.19j) in the middle to arrive at 
dm with m = 2. 
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In term (C) we use £ = II • /C + A4 and we expand the resolvent in the middle once more. 
The result is 

(C) < 4( (CI) + (C2) + (C3) + (C4) 

with 

/ \* 1 ~ V 

(CI) := (...) ^ ^ 

V J V 2 + P V 2 + P 

/ \* 1 ~ 1 ~ z> 
(C2) := (...) — -(pf + fp)^ M- ' 



P 2 + P V 2 + P V 2 + P 



(C3) := ( ■■■ ) ^ — n-/c- p 



2)2 + p V 2 + P 



P 2 + P V 2 + P V 2 + P 

Term (CI) is explicit from f)l 1.23|) after estimating one of the resolvents in the middle by 
P~ L < 1. Term (C2) is split into two terms, 



(C2) < 2(JC21) + (C22) 
with 

/ \* 1 1 ~ P 



V 2 + P V 2 + P V 2 + P 



V 2 + P V 2 + P V 2 + P 

In (C21) we use fTTTSjW uin 1 finally fTT^jl with VV = I, m = k = 2. In (C22) we use (frTTSjl 
first then (111.201) with W = P, m = k = 2. In the term (C3) we use (111.201) then (111.2411 with 
m = 1 together with the estimates ()10.14j) used for A4. 
Finally, for the term (C4) we estimate 

(C4) < 2((C41) + (C42) 

/ \* 1 1 ~ P 

C41) := ■■■) -—Ve^—U.K- 



V 2 + P V 2 + P V 2 + P 
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(C42) :=(■■•) —^—£V^ n ■ K- D 



V 2 + P V 2 + P V 2 + P 
In (C41) we first use ()11.19j) to arrive at 

(C41) < c( ■■■)*G 2 ^ U V . (11.29) 

V J V 2 + P V 2 + P 

For the term (C42) we first observe the following inequality 

Lemma 11.4. With the notations above 

1 



£V < bO(G*) . (11.30) 



V 2 + P 

Proof of Lemma \H~4\ We can write V — II • N + M, with Af,M = 0(G 2 ). Therefore 

3 3 

sv = (n • k + M)(n ■ + M) = n^-rur 1/2 e>(G 4 ) + u i°( G *) + b 1/2 o{G 4 ) 

j,k=l j=l 

after commuting II through using f)l 1 . 13[) and (jl0.14j) . Therefore f)l 1.30|) follows from ()11.17|) 
and (ITTTTBT) . □ 

Armed with Lemma fl 1.41 we see that 

(C42) < (C43) :=cb( ■■■)*G A ^ — h }C ~ V . (11.31) 
V / V 2 + P V 2 + P 

Since G 2 < G A < G s and b > 1, it is sufficient to estimate (C43) that will complete the 
estimate of (C41) from (jll.29j) as well. 

To estimate (C43), we first separate the term II • K, into terms containing H3/C3 and 
n_i_/C_i_ := rix/Cx + II2/C2 by Schwarz' inequality. We arrive at 

1 = „,/,.^ , V , / \* A 1 ~ ~ V 



(C43) < c( ■ ■ ■ Vg 4 ^ — n 3 (6 1 / 2 /C 3 ) ~ + c&( • ■ ■ Vg 4 ^ — njc 



£> 2 + P T? 2 + P v / p2 + P £>2 + P 

The first term is explicit from (|11.2(i|) with >V = n 3 , m = 2, k = 1, using that & 1/2 ^ 3 = 0(C). 
In the second term we estimate G 8 < b A in the middle, then estimate one of the resolvents 
by P' 1 < 1 and first use (111.201) and finally (lll.24jl with m = 5, together with (110.1411 for K. 
This completes the proof of (jll.6|) . 

Now we prove (jll.7|) . We need the following lemma. 
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Lemma 11.5. With the notations above, we have 



tr (•••) -(«,«) < ce^, (11.32) 

tr[---)<p£^ M ~ V (u,u) < ce- cV ~ b . (11.33) 



V 2 + P V 2 + P 

Proof. For the proof of both inequalities (jll.32|) we write £ = K, ■ 11 + Ai, then we separate the 
terms by a Schwarz' inequality and we use (jll.27|) and (jll.28|) . respectively, with appropriately 
chosen U. □ 

For the operator on the left hand side of (jll.7|) we use a resolvent expansion and a Schwarz' 
inequality to obtain 

{■■■yv^Tp<z{(D) + (E) + (F)) 



with 



(D) := (...)V>— L_ 
V J ^ V 2 + P 

f \* V 1 



V 2 + P V 2 + P 



V 2 + P v 2 + p' 

The estimate of (D) is trivial by T> = T> + £ applying a Schwarz' inequality and using (jll.27|) 
and (J11.32|) for these two terms, respectively. 
In term (E) we use 

V 1 

— V = I - P— (11.34) 

V 2 + P V 2 + P K J 

and separate it by a Schwarz' inequality: 

1 9 / \* 1 _ 1 



(E)<2 ... ^___ + 2P 2 ... £ 



V 2 + P \ J V 2 + P x> 2 + P 

For the first term we can use (jll.32j) . for the second one we use <p < 1, ()11.18j) then (jll.2z 
with m = 2. 

Finally for the term (F) we write 



(F) < 2((F1) + (F2) 
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with _ 

(Fi) :=(■■■) <p-=—iL.K- 



V 2 + P V 2 + P 



V 2 + P v 2 + p' 

For (Fl) we first estimate (p < 1, then use f)11.16|) and (jll.24j) with m = 1 together with 

fnrrm 

For (F2) we need one more resolvent expansion: 



(P2) < 2^(F21) + (F22) 



with 

1 . . V 



V 2 + P V 2 + P 
V ~, 1 . V 



(F22) := (..X(p— (VS + SV)^ M- 

v ' \ ) Y V 2 + P y J T>2 4- P ' 



£> 2 + P v P 2 + P V 2 + P' 

We split (F21) further by using V = V + S: 

f \* V V f \* 1 X> 
(F21) < 2( • • • ) y?^ .M^ + 2( • • • ) tpE^ M- 



V 2 + P V 2 + P V / V 2 + P V 2 + P 

The first term is explicit by (jll.23j) after ip < 1 and estimating T> 2 by the resolvent. The 
second term was estimated in f)11.33j) . 

Finally, to estimate (F22), we use again ()11.34|) and we split it as follows 



(F22) < 3 (...)' V £' M 



V 



V 2 + P V 2 + P 



V 2 + P x> 2 + P £> 2 + P 

P P V 
+3 ■■■ ip— — -8^ M- 



V 2 + P V 2 + P V 2 + P ' 
The first term is estimated in ()11.33|) . For the second term we use <p < 1 and ()11.18|) first, 
then (jll.26j) with W = 1, m = k — 2. For the third term we use y? < 1 and first, then 

UU} with VV = P, m = fc = 2. 

This completes the proof of (jll.7j) . □ 
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A Proof of the technical lemmas 



A.l Proof of Proposition W7I\ on the tempered lengthscale 

Proof. We recall the definitions of Bi,(x),bL(x) from (j2.1j) . (|2.2j) and we notice that Bl{x) 
is increasing in L, while bi{x) is decreasing. Since B and its derivatives are locally bounded 
and B is not constant everywhere, we easily obtain that the sets appearing on the right hand 
sides of (|2.3|) and (|2.4|) are non-empty and bounded. Therefore L m and L c are positive finite 
valued functions. 

We notice that if Bu x ){x) > L~ 2 (x) then L m (x) < L(x) < L c (x), i.e. L c (x) = L v (x). We 
claim that L c (x) = L v (x) implies that (|4.2|) and ()4.3|) hold even if L(x) is replaced with L c (x) 
which is a stronger statement. The validity of this stronger form of ()4.2j) follows directly from 
(J2.4j) . This also implies that 

B Lc{x) (x) - b Lc{x) (x) < 2L c (x) ■ sup{|VB(y)| : \x - y\ < L c (x)} < 2b Lc{x) {x) , 

therefore &l c ( x )(^) > %Bl c (x){ x )> i n particular B(y) 7^ and n(y) is well defined for all y 
with \y — x\ < L c (x). Thus f!4.3|) with L(x) replaced with L c (x) follows from (|2.4|) . We also 
proved that if L c (x) = L v (x) then the suprema in (|2.3|) and (|2.4jl are actually maxima by the 
continuity of B. 

Finally, we have to show that L(x) is tempered. Notice that it is sufficient to show that 

]x - yl < m ^i<m (A ,) 

for any x, y e R 3 because the inequality L(y)/L(x) < 2 easily follows from this. To see it, we 
assume that L(y) > 2L(x). Then \x — y\ < L(x) implies \x — y\ < L(y), so using (jA.l|) with 
x, y interchanged we arrive at a contradiction. 

Now we show that (jA.l|) holds. Let x, y be two points with \x — y\ < L(x) = ~L c (x) and 
we have to show that L(x) < 2L(y) = L c (y). This is obvious if Bu x ){y) < L(x)~ 2 , since then 
L(x) < L m {y) and L m (y) < L c {y) by definition. 

Thus we can assume that B L ^(y) > L(x)~ 2 . Since \x — y\ < L(x), we know that 

{z : \y-z\< C {z : \x - z\ < 2L(x) = L c (x)| , (A.2) 

thus B L ( x) (y) < B Lc ( x) (x), hence B Lc{x) (x) > L(x)~ 2 > L c (x)~ 2 , i.e. L c (x) > L m (x), so 
L c (x) = L v (x). 

We will now check that for 7 = 1, ... 4 

L(z;p sup {|v^|B(z)| J : \y-z\<L(x)}<b L{x) (y) (A.3) 
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and 



L(x) 7 sup | V 7 n(z) 



\y 



hold, which will imply L(x) < L v (y), hence L(x) < L c (y). But as we showed above, L c (x) = 
L v (x) implies that (J4.2j) and (|4.3|) hold with L(x) replaced by L c (x). From (|A.2|) and L(x) < 
L c (x) we therefore immediately conclude (|A.3|) . (|A.4|) . □ 



< L(x) [ < 1 



(A.4) 



A. 2 Proof of the covering Lemma 15.41 

Introduce the notation D* x := B (x, 40£(x) j and D* := D*.. Let S be any compact subset of 

R 3 . First we show how to find a finite set of points within S so that the balls Di cover S and 
they enjoy the finite overlapping property. Let D x := B(x,£(x)/20) and we cover S by the 
collection of balls D x , x G S. By compactness, we can choose points {x a } C S, with a finite 
index set a G A, such that the balls {D a } aG A cover S. Now we discard certain points from 
the collection {x a } and relabel the rest by {xj}. 

Let X\ be the point with the biggest value £{x,\) among all values {£(x a ) : a G A}. 
Then let x 2 be the point with the biggest value £(x 2 ) among all values £(x a ) such that x a G 
R 3 \ D\. Then let x 3 be the point with biggest value £(x 3 ) among all values £(x a ) such that 
x a G R 3 \ (Di U -D 2 ), etc. until all covered by D^s. This selects a subcollection of 
the points {x a } and they are relabelled to Xi, x%, 

We claim that the collection of D^s cover S. Consider any y G S, then y G D a for some a. 
But x Q is covered by some D^. We choose the smallest such index i. By the maximality of the 
radii in the selection procedure, we know that £(x a ) < £{xj) , so \y — Xi\ < \y — x a \ + \x a — Xi \ < 
{£{x a ) +£{x l ))/20 < £{xi)/10, hence y G D f . 

We claim that the union of the D* balls have the finite covering property with a sufficiently 
big universal N. From construction, the balls Df := B(xi,£(xi)/40) are disjoint. Fix a point 
y G R 3 and let / be the set of indices i such that y G D*, i G /. Choosing e < 1/40, we see 
that y G -D*, i.e., \xi — y\< 40£(xj) implies 1/2 < £(y)/£(xi) < 2 for all z G J. Hence the balls 
z G /, all have radius at least £(y)/80 and they are within a ball of radius 81£(y) about y. 
From their disjointness it follows that their number is universally bounded, i.e. the number 
of -D*'s covering any y is bounded by a universal number N. 

This completes the construction of the covering balls for any compact set S satisfying (i) 
with a universal covering property. 

We denote by P(S) the points {x\, x 2 , ■ ■ .} obtained in this procedure and note that P(S) C 

S. Let H(S) := U ie p ( 5) A and H*(S) := U i6 p( 5 ) A* 

Now we show how to choose points in the whole space. Fix an arbitrary point x, and let 
A k := {y : 4 k £(x) < \y — x\ < 4 k+1 £(x)}, k = 1, 2, . . ., be a sequence of annuli. Clearly 
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D x U \J k A k = R 3 . For each annulus we construct the points P(A k ) defined above and we let 

oo 

P:=P(D x )u[jP(A k ). 

k=l 

This will be our final set of points {x\,X2, ■ ■ ■} after relabelling. It is clear that the balls 
Di = B(xi,£(xi)/ 10), Xi G P, cover the space. 

Next we prove the finite covering property for the balls {D Xi : X{ G P}, i.e. that the 
number of balls that cover any given point of R 3 is universally bounded. We need a lemma 
whose proof is given later. 

Lemma A.l. Fix any point x G R 3 . 

(i) Let L k : = sup{£(w) : u G A k }, k>\, then L k < A k+1 £(x). 

(ii) If y G A k , then y H*(A m ) for any \m — k\ > 5. 

Recall that for each m every point in R 3 is covered by at most N balls D with center 
z G P(A m ) and similarly for balls with center in P(D X ). Hence (ii) of Lemma fA. II shows that 
any y is covered by at most 12N balls with center from P. This completes the proof of the 
finite covering property of the balls {D Xi : ij 6 P}. 

Finally we show property (ii) of the Definition ^. HI If DidDj ^ 0, then £j, Ej are comparable 
by (|4.1|) . Therefore D* covers Xj, but any point is covered only by finitely many D*'s, hence 
Dj can be intersected by finitely many D^s. 

Proof of Lemma \A.ll (i) Suppose that there exists u G A k with £(u) > 4 k+1 £(x). Then 
\x — u\ < A k+1 £(x) < £(u), hence £(u) < 2£(x) by (|4.1j) which is a contradiction. 

(ii) Suppose that there is a point z G P(A m ) such that y G D*, i.e., \y — z\ < 40£(z). Using 
IO this implies £{z) < 2£{y) assuming e < 1/40. Hence \y-z\ < 80£(y) < 80-4 fe+1 £(x) by (i). 
But z G A m , so \y—z\ > (4 m — 4 k+1 )£(x) which is a contradiction if m > k+5. Suppose now that 
m < fc-5. Then \y-z\ > (4 fc -A m+1 )£(x) which contradicts to \y-z\ < A0£(z) < A0-A m+1 £(x). 
□ 

A. 3 Proof of the localization Proposition 17.11 

As a preparation for the proof we define a distance function on the collection {xi} iG j obtained 
in Lemma [5.41 Let i G / be a fixed index. We define the following compact sets successively 

S (i) := Di , 
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Sfc+iOO : = IJ D j , 

and we denote 

m k := card{j El : Dj Pi Sfc(z) 7^ 0} . 

Lemma A. 2. (i) The sets So(i),Si(i), . . . are increasing. 

(ii) Let u E Skii), v E~ int(Sk+i(i)) , then \u — v\ > M{u). 

(iii) \J k 3 k (i)=R 3 . 

(iv) rrik < N k+1 with the universal constant N from Lemma \b.J\ 

(v) For sufficiently small e and for any nonnegative function G such that G(x) and G(y) 
are comparable whenever \x — y\ < s~ 1 £(x), we have 

sup G <2 k sup G . (A.5) 

S k (i) S (i) 

Proof. For simplicity, we omit i from the arguments since i is fixed. 

(i) Since the balls {Dj} cover R 3 and Dj C Dj, we see that Sk C Sk+i- 

(ii) Let u E Dj for some j. From \u — Xj\ < £{xj), it follows that £{xj), £{u) are comparable. 
Since Dj C 5*^+1, we have \v — Xj\ > 10£(xj), so \v — u\ > 9£(xj) > A£{u). 

(iii) Suppose that 5* := |J fc 5*^ is not the whole R 3 and select a point z E OS. Then we can 
find a sequence of points Zk E dSk converging to z such that \zk — z\ monotonically decreasing 
(for example, we can choose the point z k E Sk closest to z). Since £{z) > 0, we see that 
\ z — z n \ < £(z) for some n, hence £(z n ) and £{z) are comparable. We have 

\z n +i ~ z n \ < \z n+l - z\ + \z- z n \ < 2£(z) < A£(z n ) 

which contradicts (ii). 

(iv) It is clear that mo < N by Definition 15. 31 ( ii) . By induction we show that m>k+i < Nrrik- 
This again follows from Definition l5.3l fii). since each Dj in the definition of Sk+i may intersect 
at most N balls from the collection {Dj}j^i. 

(v) Straightforward by induction on k and by the definition of Sk(i). □ 

This lemma gives an integer valued distance on the collection {xj}j e j: 

dij := mm{k : Xj E Sjfe(i)} • 

Clearly da = 0, and dy + djk > dik, but the distance function is not symmetric. However, we 
have 
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Lemma A. 3. For sufficiently small e the distance function satisfies 

dji < 7dij + 1 . 



(A.6) 



Proof. The proof goes by induction on the value of dij. If d^ = 0, Xj G -Dj, then Xi G Dj, i.e., 
dji < 1. Suppose that ()A.6|) is proven for all pairs with d^ < d and let now = d + 1. 
Then there exists m such that x m G Sd{i), di m < d and Xj G -D a for some index a with 
\x a ~ x m \ < £ a + I0£ m . If e is sufficiently small, the radii £ a , Ej and £ m are comparable, and it 
easily follows that dj a < 3 and d am < 3. Therefore dp < dj a + d am + d mi < 6 + 7d+l < 7^ + 1. 
□ 



For every i G /, k G N we define 



u 



j : Z»jnS fc (i)^0 



(A.7) 



Notice that supp(w[ ) C Sjfe + i(i), and «j/ = 1 on Sk(i). Moreover 



(A.8) 



To see this, we notice that every x is covered by not more than N balls Dj and only these 
support those fy's which do not vanish at x. Moreover £j is comparable to £{x) for all these j 
indices, hence HV^Hoo < c£{x)~ x . 

In the rest of this section we set 

R f = R(f) := (T + /)~\ Ri\f\ := (T t + f)' 1 

for simplicity, in accordance with the notations ()6.2|) . 
Proof of Proposition \7.1\ We start with an auxiliary lemma. 

Lemma A. 4. For any number \x > and rea/ function \ on R 3 , 



1/2 



1/2 



< 1 



< c 



p- 1/2 |v x | 



(A.9) 
(A.10) 



Proof of Lemma \A.4 The first inequality is trivial by inserting P + fi<T-\-P + fi. For the 
second inequality we use 

[T,x]=A*[A,x} + [A*,x}A, (All) 
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and it is sufficient to estimate one of these terms; 



< 



p 1 / 2 A* 



[A, X ]R 



1/2 



Using ||M|| = || MM* II 1 / 2 , we obtain that the first factor is bounded by 1 (again, using 
A* A = T<T + P + fi). For the second factor we need pointwise commutator bounds 



[A, x }[A,xT < c |Vx| 2 , [A*, X }[A*, X ]* < c |V X p 



(A.12) 



that follows from ()7.2|) . 

Hence, we estimate the second factor as 



< c 



1/2 



< c 



R l/2 poV2 



1/2 



P~ L \Vx\' 



1/2 



and we use flA.9|) . □ 

The key lemma is the following (recall the definition of Xi from Section l5~2j) : 
Lemma A. 5. For sufficiently small e we have the following estimates for any i,j(zl 

XiRp+pOjFHjRp+pXi < c{Ac e) 2{d ^ + Ff{P l + ^XiRp^Xi , 
XiRp +l tF 2 Rp +ll Xi < cF 2 (P t + ii)- l XiRp + nXi , 



XiRp+aXi < (Pi + v) 



-2 



XiAR 2 P+u A* Xi < c(P % + ^Y 1 . 



(A.13) 
(A.14) 
(A.15) 
(A.16) 



Proof of Lemma \A.5[ For brevity, we denote R := Rp+^. First we show f)A.13|) . We assume 
that Eq is small enough so that F(xA and F(xj) are comparable as long as Dj D Di ^ (see 

(EH))- 

We first consider the case dij < 1. Then F 2 < cF 2 < cF 2 on the support of 9j. We also 



use 



2 < c(Pj + ii)-\P + n)< c(P 4 + n)-\P + fi) . (A.17) 
Hence (IA.13J) follows from 

XiRejF^jRxi < cF 2 XiRe 2 R% < cF 2 {P, + fi)~%R(P + ^)R Xi 

and finally we use (|A.9J) to estimate R(P + fi)R < R. 
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To prove (jA.13|) for d = dy > 2, we recall the definition of the functions ujj. (|A.7|) . For 
brevity, we omit the superscript i. We successively insert the functions Ui,u 2 , ■ ■ • where 
d = d. 



i.r 



XiROj = XiUiROj = XiR[T, ui]R9j = XiR[T, Ui]u 2 R9j = XiR[T, Ux]R[T, u 2 ]R6j 
= ... = XiR[T, Ul ]R[T, u 2 ]R. . . R[T, u d ^]R9 3 . 

We used that u\ = 1 on the support of Xh u k+i = 1 on the support of Vw/c and suppf/u^-i) n 
supp(0A = 0. Therefore we can first estimate 6jF 2 6j < cFj9 2 , then use the succesive insertions 
to obtain 



XiR6jF 2 9jRxi 
< cF?XiR 1/2 



(A.18) 



d-l 



n (r^t^r 1 



/2 



k=l 



d-l 



/2 



L fc=l 



R 1 ' 2 



Xi 



First we use that 

R l/2 6 2 R l/2 < c (p. + M )-l jR V2(p + < c (p. + 

by (jA.17j) and (|A.9|) . Then we use (|A.10|) to estimate the commutator norms and we use (jA.8|) 

to get 

P- 1/2 \Vu k \ <ce 5 / 2 <e 

oo 

for sufficiently small e. We obtain 

XiR9jF 2 9jRxi < c(c Q e) 2 ^^F 2 (P 3 + ^XiRXi (A.19) 

By ()A.5|) . we see that P/(-Pj + At) -1 < lG^F^Pj + /i)" 1 because part (v) of Lemma lA~2l applies 
both to the function G = F and G = (P + n) 1 ^ 2 . This completes the proof of (|A.13I) . 

To prove (jA.14|) we insert a partition of unity 

XiRF 2 Rxi = Y.XiRdjF^jRxi . 



We use (lA~T3j) . (iv) of Lemma lA~2l and that 



^(4c e) 2(dl ^ 1)+ < 1 + N + ^(4c e) p A^ p < N + 2 

jG/ p=l 



(A.20) 
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if e < Eq, where the universal constant N is from Lemma [5.41 This proves (|A.14j) . 



The proof of ()A.15|) is straight-forward by applying (jA.14j) with F = 1, 



(A.21) 



and then using i?p +M < (P + fi) 1 which is bounded by c(Pj + fj) 1 on the support of x% since 
P and Pi are comparable on this set. 

For the proof of (jA.16|) we insert x% that is identically 1 on the support of Xi an d use 



We can remove Xi ari d use ARp +fl A* < 1 to finish the proof. □ 

The next lemma is a strengthening of (jA.13|) in Lemma [A. 51 Notice that in (jA.13|) we lost 
a resolvent, and the right hand side is not locally trace class in the high momentum regime. 
The following lemma remedies this: 

Lemma A. 6. For sufficiently small e 



Proof of Lemma \A.(A For simplicity, we let R := -Rp+ M and R{ := Ri[Pi] in this proof. We 
start with the identity 



(EH) 



XiAR 2 P+fM A*Xi = Xi A XiRp +tl Xi A *Xi < c{Pi + n) 1 XiAxiRp+ i iXiA*Xi ■ 



d]R P ^d 3 F 2 d 3 R P+ii d] < c^ef^+Fffl (i2»[P 4 ] + P^R^A^AR^O 2 . (A.22) 



%R = RiXi + Ri (xi{Pi -P-») + A*[A, Xi} + [A*, %}A)R , 



(A.23) 



since A and Aj coincide on the support of % D Y Q7.3J1 . After a Schwarz' inequality 



XiR9jF 2 9jRxi 

< jRaiOjF^XiRi + RiXi(Pi -P- ^)R9 j F%R{P i -P- n)%Ri 



(A.24) 



+RiA*[A, x i }R6 3 F 2 9 j R{A, XifAiRi + Ri[A*, Xi\AR9jF 2 9jRA*[A*, Xi}*Ri 




The first term is estimated as 



RiXiOjF^jXiRi < cF 2 R 2 l{d i:j < 1) , 



(A.25) 



since F < cF { on the support of 9jXi- 
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Since Xi = 1 on the support of Xi, we can freely insert Xi m the last three terms of (jA.24|) 
replacing R9jF 2 9jR with XiR@jF 2 @jRXi everywhere and apply (jA.13|) after multiplying it by 
Xi from both sides: 

\,/^/^ ; //\, < c^ef^-^F^P, + fi^XiRXi ■ 

In the second term of (jA.24|) we use R < (P + fi)^ 1 and that P + ji and Pi + \i are 
comparable on the support of %. Hence 

RiXi{Pi -P- ll )nO J l--0 J h'(l', - P - ti)xiRi < c(Ac s) 2 ^- 1)+ F 2 R 2 , (A.26) 

using (Pi — P — n) 2 < (Pi + /i) 2 on the support of x%. 

In the third term of (jA.24|) we again estimate R < (P + fi)^ 1 < cP^ 1 on the support of 
Xi, we use (jA.12|) and that | [A, Xi][A, Xi]*\ = |Vxi| 2 < ce 2 Pi^ 2 to obtain 

RiA*[A, XijROjF^jRlA, Xi}*AiRi < c(4c e) 2d - Ff Pr l R^^AiRi . (A.27) 

Finally, the fourth term of (jA.24|) satisfies 

Ri[A*, Xi] AR6j F 2 6j RA* [A* , %}*R % < c^e) 2 ^- 1 ^ F^P^R^A*, Z]A Xi RXiA*[A\ Z]*R % 

< c(Ac e) 2d ^F 2 R 2 , (A.28) 

since x% = 1 on the support of Xi, we can omit it, and we used ARA* < 1 and (jA.12|) . Lemma 
ESJ follows from (|A3^ - (jA~28j) using that Bcli = Oi- □ 

Finally, we complete the proof of Proposition 17. II We use R = Rp+^ for brevity. We insert 
three partitions of unity and perform a weighted Schwarz' inequality 



RF 2 R = 2 iR0jF 2 6 j Re 2 k (A.29) 

(e^- 1)+ -^- 1)+ e 2 Re j F 2 e j R6 2 + e^- 1) +- {dk 3~ 1) +6 2 k R6 j F 2 8 j R8 2 k 

4 L~c T \ S 



i,j,kel 

< 



i,j,k&I 

Using (jA.6|) we see as in (jA.20|) that 

^ e Ww-l)+ < y^ e (dj*-8)+/7 < ^9 
k&I k&I 

if e is small enough. Here N is from Lemma [5.41 Therefore (|A.29J) implies 

RF 2 R < cN 9 J2 £ ~ {dlJ ~ 1)+e i Re j F2d j Re i ■ 
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We use (jA.22j) and sum up the index j similarly to (|A.20J) with a possible smaller e 
RF 2 R < cN 8 J2( 16c l £ ) idij ~ 1)+F i i ( R i + P t ~ 1 RiA*(p 2 i ARi)0* 
< cN 10 FM (r 2 + P-'RiA^ARi) Oj . □ 

i 

A. 4 Proof of Lemma 18.21 on the magnetic coordinates 

We give the construction of the new coordinates and conformal metric but we do not follow 
the explicit bounds along the proof as they easily follow by scaling. 

We consider z G V fixed in the proof and omit the notation z in the sub- and superscripts. 
We define the function 

«(r):=-~tog|B(^(r))| (A.30) 

where k(t) G C 3 (R). At each point <f(r) we consider a small spherical cap of the sphere S(t), 
going through <f(r), orthogonal to the field line if and having curvature |/c(t)|. The different 
spheres should curve in a direction determined by the sign of k(t): positive curvature means 
a sphere with outward normal pointing in the direction of if. We consider now a small 
cylindrical tubular neighborhood, Af := {x : inf r dist(x, fij)) < 10£}, along the C 5 -curves 
ip, in which the spherical caps define a C 3 -foliation. Since B is extended (D, i^)-regular, such 
neighborhood exists if e < e(K) is sufficiently small. The foliation naturally extends the 
function r onto Af such that it is constant on the leaves and r G C 3 (R 3 ). 

This foliation can be extended to the whole of R 3 in the following way. Since B = 
outside of D = B(zo, £), for some zo G R 3 , in case of \z — zq\ > 2£ we have k(t) = and if is 
a straight line, so the foliating spherical caps are parallel flat discs and they can be trivially 
extended to a foliation of R 3 with parallel planes. 

Now we consider the case | ^ — ^ 1 < 2£, where we only know k{t) = for \t\ > 3£. In 
this region the spherical caps are again parallel flat discs and they can be extended to parallel 
planes. That leaves a parallel slab unfoliated between the planes passing through <p(—3£) and 
<p(3i). The width of the slab is (6 ± cKe)l 

We consider the smooth function F : R 3 h> R of the form F(x) = A(noo • x), where 
A : R — > R is smooth and is chosen such that F(ip(r)) = r for |r| > 3£, ||A' — l||oo < cKel 
and || A 1 ' 7 -' \ oa < cKei - " 1 , 7 = 2, 3, 4. The level sets of F define a parallel foliation of R 3 which 
coincides with the previous foliation outside of the slab. 

Let x be a smooth cutoff function supported on the ball D := B(zq,6£), || V 7 x||oo < c£~ 7 , 
7 = 1, . . . , 4, and \ = 1 on B(z , 5£). Since \z — Zq\ < 2£, we note that D C Af, so r is already 
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defined on supp(x). We define the function 

t : = X t + (1 - X )F . 

An easy calculation shows that t G C 3 (R 3 ), ||Vi — n^H < cKe and the level sets of t define a 
regular foliation of R 3 . This is clearly an extension of the foliation given by r on D and the 
leaves are planes on D c . Moreover, if we define a smaller tubular neighborhood M of the central 
field line as Af := {x : inf T dist(x, <p(r)) < 2i}, then we note that Af C B(z , U) U {|r| > 3£}, 
therefore t = r on Af. Let iV := ||V£|| -1 V£ be the unit C 2 -vectorfield orthogonal to the 
foliation. We remark that the integral curves of iV typically do not coincide with the field 
lines except on the field line ip and in the region far away from zq. 

Armed with this foliation, we introduce new coordinates on R 3 . On the plane V we choose 
Euclidean orthonormal coordinates £i,£2 with origin at z and dual to the basis {pi,P2}> i-e. 
x — z = £ipi + (,2P2- Clearly (d^,, <%.) = 5j k for j, k — 1,2. For simplicity we set dj := 8%. and 
Vj := V dj in this proof. 

We extend the coordinate system £1, £ 2 defined on the plane V by setting £1, £ 2 constant on 
the integral curves of N. It is easy to check that £1, £2 £ C 3 (R 3 ). Together with t they define 
a regular set of coordinates on R 3 . The central line is given by (0,0, t) in these coordinates. 
Let b(t) := |B(0,0, t)\ G C 3 (R) be the strength of the magnetic field along the central line, 
note that b(t) is comparable with b for all t by (j5.1|) . We define £3 = £3(2) to be the solution of 
^£3 = {bfy/b} 1 / 2 with £3(0) = 0, clearly £ 3 £ C 4 (R). We reparametrize the coordinate t with 
£3 := £s(t). In this way we defined a new coordinate system, {£i,£2,£s}, with origin at z. We 
shall view the coordinates £i,£2,£3 as C 3 functions of x G R 3 and whenever the dependence 
on z is relevant, we use the notation £ 2 = (£i,£2,£|)- It is easy to check that the function 
(z,£) 1 — > £ 2 belongs to C 3 (R 5 ) and the vectorfields di = % are C 2 . 

We note that in the trivial case, \z — zq\ > 2£, we simply have £ 2 (x) = i?'(a; — z), where 
R := [pi|p2| n oo] is the 3 by 3 matrix with columns pi,p2,n OQ , and all statements of Lemma 
18.21 are trivial with Q = h = 1. 

From now on we shall assume that \z — z \ < 2£. The relations (|8.3j) - (j8.4|) and (|8.6j) . i.e., 
the fact that ds 2 has no d£j d£3 (j = 1, 2) components follow directly from the construction. 
From the regularity of the magnetic field (Definition 15. ip it easily follows that the Jacobian 
of the change of coordinates x 1— > £ z (x) is close to the matrix R l and it varies regularly in z. 
This proves that the function (z, £) 1— > x 2 (£) is well defined and C 3 (R 5 ), it also proves (|8.5|) 
and (J8.7|) by the inverse function theorem. 

The metric is diagonal in the £ coordinate system on the plane V, i.e. for t = 0. The key 
point is to show that it remains diagonal within the tubular neighborhood Af. The diagonal 
metric elements will define the functions Q and h. 

We derive a differential equation for the metric components '■= (dj,dk) g G C 2 (R 3 ) 
where j,k = 1,2 within N. We have d t g jk = d t {dj,d k ) g = (V t dj,d k ) g + (dj, V t d k ) g . Using 
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that d t ,dj,dk, are coordinate fields, i.e, have vanishing Lie derivatives we have Vtdj = Vjd t . 

— 1/2 

Recall that N = g tt dt is the unit normal to the spherical foliation, where gu '■= (dt,dt) g . 
Then V jd t = g]( 2 V jN + dj{g]( 2 )N and therefore we have 

dtQjk = ^/ 2 [(V,iV, d k ) a + (9,-, V k N) g ] = 2g\l 2 K 3k , 

where Kj k is the second fundamental form of the leaves of the foliation. For a sphere immersed 
in R 3 we have Kj k = ngj k , where k is the curvature. We recall the choice of k from (jA.30|) 
and that t = r on Af. Thus 

dtgjk = 2gl/ 2 n(t)g jk . (A.31) 

This proves that since g±2 is zero on the supporting plane t — 0, it is zero everywhere 
in M. It also proves that gu = g 2 2 everywhere in M since they satisfy the same equation 
and initial condition. The same relations trivially hold for the region > KM, where 
<?n — #22 — 9\2 — 1- Moreover, we define Vt := g^^ 2 G C 2 (R 3 ) and we obtain that within J\f 
as well as in the regime > 10£ the conformal metric can be written as 

dsl = d^ 2 + d£ 2 + n 2 gtt dt 2 = d^ 2 + d£ 2 + n 2 gtt f(&y 2 de 3 

_____ 1 /2 

using the definition of / and £ 3 . This proves (|8.8j) with h := f2p t / /fe). Since the new 
coordinates form an orthonormal system for |^| > 3£ and also for > 10£ modulo a change 
of variables in the third direction, the identities ()8.9|) and (|8.12|) . respectively, follow from the 
definitions. 

Along the central line we have g tt (ip(t)) = 1. Thus (jA.31|) . ()A.30|) and gu = 1 for |r] > 3£ 
implies that gu = b/\H\, i.e., Q = f{C,s) and h = 1 along the central line. Then (J8.5|) implies 
()8.10|) . The global bounds (j8.13|) - (j8.14|) also follow from the smoothness of the contstruction, 
i.e. from ()8.5|) . The details are left to the reader. 

Finally, the orthonormal basis {ei,e2,es} in ds^ is defined by first constructing e[ := 8%, 
e' 2 := 82, e' 3 := h~ x 8^ which are automatically orthonormal apart from the region {£ : |£ < 
< 9£, 1^1 < 3£}. On this region we apply a Gram-Schmidt orthonormalization procedure 
to obtain {e\, 62, 63} from {e[, e' 2 , e' 3 }. □ 



A. 5 Comparison of operators on equivalent L 2 -spaces 

Let dfj, and dv be two positive measures on R rf , that are mutually and uniformly absolutely 
continuous, i.e. dv(x) = F(x)dfi(x) with a a positive bounded function F with bounded 
inverse We let 

r — HFii iip-iii - maxF 

L-F •— \\-r do F 00 — : ■ 

mm r 
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Consider the spaces L 2 (d/i) and L 2 (dv) = L 2 {Fdfi) and let A be any operator defined on 
L 2 (d/i). Since these two spaces are the same as sets, we can consider A acting on L 2 (Fd/i) 
as well. We denote this operator by Ap. Let (•, •) := (-, -)-L 2 (d^) and and (■, -) F '■= (•, -)l 2 {faij)- 
Similar convention is used for || • || and || • \\p and for the traces over these L 2 -spaces: Tr := 
Tr L 2 (dM) and Tr F := Tr L 2 {Fd ^ 

Lemma A. 7. Let A be a Hilbert- Schmidt operator on L 2 (d/z) with a kernel A(x,y). Then Ap 
is also Hilbert- Schmidt on L 2 (dv), and the kernels of these operators satisfy 

A(x,y)=A F (x,y)F(y) . (A.32) 

Moreover, the diagonal kernels of A* A and A* F Ap are comparable: 

C F 1 \\F\\^(A*A)(x,x) < (A m F A F y x ,x) < CfWF-^A* A){x,x) . (A.33) 

Furthermore, if < a < A* A < (3 for some constants a, (3, then 

C F l a< A* F A F < C F (3 . (A. 34) 

If, in addition, A is of trace class, then so is Ap and their diagonal kernels satisfy 

A(x,x) = A F (x,x)F(x) . (A.35) 

Proof. Recalling the conventions at the end of Section the identities f|A.32|) and (|A.35|) are 
obvious. For (jA.33|) we estimate 



(A F A F )(x,x) = J A F (x,y)A F (y,x)F(y)dfi(y) 

= F~ 2 (x) J \A(x,y)\ 2 F(y)d^y) 

< C F \\F- l \\ oc (A*A)(x,x) . 

The lower bound is proven similarly. 
For ()A.34|) we notice that 

(i/),A* f Af*I>) f = \\F l ' 2 A^\\ 2 < (maxF)II^H 2 </?(maxF)||^|| 2 

and 



< \\F~ 1/2 n 2 F < (maxF- 1 ),,^,^ , 
which proves the upper bound. The proof of the lower bound in (jA.34|) is similar. □ 
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Lemma A. 8. Let be a finite collection of closed operators on L 2 (dfi), let W\, W 2 be 
nonnegative functions on H d . Then 

Tr F ^J2( A kT F (^k)F + W 1 -W 2 y\ < J Tr^J2 A l A k + W 1 -C F W 2 ^j_ | . (A.36) 

k k 

Proof of Lemma \A.S\ By the variational principle 

Ti F ( K Y J ( A k) F {A k ) F + W 1 -W 2 )_ = inf {Tr F (22(Ak)* F (A k ) F + W 1 -W 2 )'Y : < 7 < l} 

k k 

where the infimum is over all finite rank density matrices 7 on L 2 (Fdfi). We can write 
7 = J2 n ^n(fn, OWn with < A n < 1 and {/„} being orthonormal in L 2 (F). 
Define the operator 7 := (minF) J2 n ^n(fn, -)fn 011 L 2 . Since 

(0,70) = (minF)^A n |(/ n ,F-V)Fr = (mm F)(F~ 1 4 > , 1 F~ 1 ( f ) ) F < ||0|| 2 , 

n 

7 is a density matrix on L 2 . Furthermore, for any A = A k 

Tt f A* f A f1 = Y, X n\\A F f n \\ 2 F = Y,U\ A fn\\l>(mmF)J2U\Afn\\ 2 

n n 

imF)J2^nTr\A*Af n )(f n \ = TrA*Aj . 



n 

(mm. 



The potential term is estimated as 

n 

> {mmF)J2*n{fn,(W 1 -C F W 2 )f n ) 

n 

= Tr(Wi - C F W 2 )^ . 

Therefore 

Tr F (^(A fc )^(A fc ) F + Wx- ^2)7 > Tr(^A*A fc + W X - C F W 2 )^ 

k k 

and (jA.36|) follows from the variational principle. □ 



67 



A. 6 Comparison of Dirac operators under a conformal transforma- 
tion 



Let Q : R 3 — > R + be a C 1 -function satisfying 

1 



(A.37) 



and 



IVOI 



< 



(A.38) 



with some constant I > 0. We define the metric ds^ := Q 2 ds 2 that is conformally equivalent 
to the Euclidean metric ds 2 . Let D be a Dirac operator in the ds 2 metric, then a Dirac 
operator in the ds 2 ^ metric is given by T>q := Q~ 2 VQ. Notice that T>q is self-adjoint on 
L 2 (ds^) £g>C 2 (see jES-lllj ). The following lemma compares certain resolvent kernels of T> and 
Vq on L 2 (ds 2 ) (g) C 2 and on L 2 (ds^) ® C 2 , respectively. 



Lemma A.9. Let P > 2 9 £~ 2 be 



a num 



her. Under the conditions \A.31 ), \A.38\) we have 



tr 



(v 2 + py 



(x,x) < 2 9 tr 



(vi + py 



(x, x) 



x e R 



(A.39) 



/ 2 



TTie /e/t hand side is the diagonal of an operator kernel on L 2 (ds 2 ) (g) C 2 ; i/ie n'<?/ji hand side 
is the diagonal of an operator kernel on L 2 ((1sq) ® C 2 . Moreover, if < (p < I is a bounded 
function then 



tr 



£> 2 + P 



V<p 2 V 



1 



£> 2 + P 



(A.40) 



L 2 



< 2 12 tr 



P 2 +P 



T>n<p 2 Vt 



m + p 



p 



(Pi + py 



[x, x) 



i 2 



for any x G R 3 . 



L 2 (ds^) ® C 2 be a unitary map given by 



Proof o/ Lemma E3 Let V : L 2 (ds 2 ) <g> C 2 
W := fi~ 3 / 2 ^. Notice that 

v = v*(n 1/2 v n n 1/2 )v , 

therefore the unitary operator Vl 1 I 2 T>qVi}I 2 on L 2 (ds^) ® C 2 is unitarily equivalent to D on 
L 2 (ds 2 ) ® C 2 . In particular, for any real function / 



f{v) = n 3/2 f(n 1/2 v n n 1/2 )n~ 



3/2 



(A.4i; 
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From (jA.41|) we obtain 

1 4 



< 



(V 2 + P) 2 ~ (V 2 + P) 2 + 3P 2 
In particular, 

1 



ft 3 / 2 



{[Q^VnQ 1 / 2 ] 2 + P) 2 + 3P 2 



(v 2 + py 



;(x,x) < 



(OV2p n OPoOV2 + py + 3P 2 



(x, x) . 



L 2 



Here the right hand side is the L 2 (ds 2 ) <S> C 2 kernel of the corresponding bounded non self- 
adjoint operator. However, the same operator can be viewed on L 2 (ds^) <g> C 2 as well, where 
it is self-adjoint. Using (jA.35|) from Lemma [A. 71 we know that the two diagonal kernels differ 
by a factor Q 3 (x). 

To conclude (jA.39|) . it is therefore sufficient to show that 
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ft- 1 / 2 



(nV2z> n nz> n fii/2 + p) 2 + 3P 2 - {pi + p) 

as self-adjoint operators on L 2 (ds 2 7 ) ® C 2 . Using 

Q^VnQVnQ 1 / 2 = Q^V^Q 1 / 2 + Q 1/2 [V n , Q}V n n 1/2 
and a Schwarz' inequality, we obtain 
(n^VntlVnn 1 / 2 + P) 2 + 3P 2 



(A.42) 



(A.43) 



> ^(Q 3 / 2 P 2 ft 1 /2)*(Q3/2p2 n l/2 ) _ 2 ^ 1 /2[p n) Q]p n Ql/2j ^[Pn, ft]^ 1 / 2 ) + p 2 



> 



ft 1/2 p^ n 3 v 2 n n 1 ' 2 - AQ}/ 2 v n [Vsi, n]*n[v n , Q\v n Q}/ 2 + 2P 2 
i 



.fiV2p4 l/2 _ 8n l/2£> n [Z> n , ^P^ft 1 / 2 + 2P 2 

Using that [Pn, fi] = ft~ 2 [P,f2]f2 and (|A.38|) we can estimate 

[VnM*[V*M <4||Vft||L<2- 7 P 



(A.44) 



so we can continue 



flA~43l) > =p 1 ' 2 



Ip4 _ Lpv 2 

8 a 16 n 



P 2 



nv* > -LnV 2 (v 2 n + pfn 1 ' 2 . 
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This completes the proof of (jA.42j) and hence (jA.39|) . 

For the proof of (jA.40j) we can use the argument above to reduce the problem to estimating 
the diagonal element of the self-adjoint operator 



T := m^Vn^'^^V^R with R : = 



1 



(fi 1 /2p f1 Ql/2)2 + p 



viewed on L 2 (ds^) (g> C 2 , where Vq is self-adjoint. The resolvent can be written as R 
Sl-^RxQ.- 1 / 2 with 



Ri 



and we can expand 



1 



V n QD Q + Pfi 



-1 



n- 1 



Rx = Q- 1 



Ri 



ivl + p 



Therefore, using a Schwarz' inequality, ()A.37|) and < (p < 1, we can estimate 
T < 2tt- 1/2 R 1 V n <p 2 V n R 1 n~ 1/2 



(A.45) 



vl + p- 



+1 



VnRxVnlV^np- 1 



v n RxP(n- 



V 2 n + P 



V 2 +P 



Qr l/2 . 



Here we used the shorthand notation (■ • -)*A for the operator A* A where A is a long expression. 

In the first term on the right hand side of (jA.45|) we use VqQ^ 1 = Vr l V^ + [Dq, Q" 1 ] and 
fTOil to obtain 



< 



vl + p- vl + p 



2 -8p 



{VI + P) 2 



and both terms explicitly appear on the right hand side of (jA.40|) . For the other two terms it 
is sufficient to show that 



RxVlRx < AP- 1 
V n RxV 2 n RxV n < 4, 



(A.46) 
(A.47) 



70 



and then the last two terms in (jA.45j) can be estimated by the second term on the right hand 
side of (EOOj) using (pPTj) . (|A~38j) . 

For the proof of (TPoT) and (TP7I) we first use < 2V n ttV Q + 2PfT 1 = 2K[ l to cancel 
one of the resolvents. The proof of (jA.46|) is then completed by estimating the other R\ by 
2P~ l . For the proof of (jA.47|) we notice that 



1 _ 2 

2 + pn- 

This completes the proof of Lemma lA. 91 □. 



VnRiVa - Vn VnnVn + pn _ 1 V Q < ^^^^ < 2 



A. 7 Proof of Lemma lll.2t apriori bound on the full resolvent 

Using (jA.34|) from Lemma [A. 71 and since the volume forms du and d/i = Q 3 dx are comparable 
at every point, it is sufficient to prove (|11.15j) - (jll.l9|) in the space L 2 (d/z) = L 2 (Q 3 dx). In 
this space V = and the components of D = = 11^ — |(divn/i, divn/ 2 , divn/ 3 ) are 
self-adjoint. We recall that 11^ was given in (j9.14j) and D was defined in general in ()9.7|) . 
Throughout the proof we will work in the space L 2 (Q 3 dx), and we adapt the notation T> = Dq, 
D = in this section. We also recall that IT, = Dj + idj with dj := |divQ ej. 

Using Lichnerowicz' formula (|9.9|) . sup^. = sup^. ||B(x)|| < cb and that all geometric 

terms are bounded by (j8.13|) and (j8.14|) . we can estimate 

V 2 > D 2 - cb . (A.48) 

We recall that £ = 1, b > e~ 2 > 1 and P = e~ 5 > 1. 

For the proof of (jll.l5j) we start with a Schwarz' inequality 

and use that \dj\ < c. We estimate one of the resolvents trivially and use (|A.48|) 

Dj(—^ Yd, < D j —^- D j < Dj— D j < Dj — t, — -Dj < b 

3 \V 2 + pJ 3 ~ 3 V 2 + Pb 3 ~ 3 B 2 -cb + Pb 3 ~ 3 D 2 + Pb/2 3 ~ 

for sufficiently small e. This completes the proof of 

The proof of (lll.lftj) is identical just we estimate (V/(V 2 + P)) 2 by (V 2 + P)" 1 . 

For the proof of ()11.17|) we first compute HjHk = DjD k + i{Djd k + D k dj) — (<9 efe c4) + djd k - 
We use a Schwarz' inequality, the estimate ()11.15|) and the boundedness of dj's together with 
their derivatives, we obtain 

u * kU j(wTpY UjUk - 2DkDj (wTp) DjDk + cb ■ (A ' 49) 
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We apply Lemma [11. II to estimate the resolvent square, using that T> 2 and D 2 differ only by 
an operator bounded by cb < Pb/2 if e is sufficiently small: 



< 



< 



Ab 2 



V 2 + PJ ~ (V 2 + Pb) 2 ~ (D 2 ) 2 + (Pb) 2 /4 



(A.50) 



We expand (D 2 ) 2 = J2j Dj + J2 j<k {D 2 Dl + B\D 2 j) and use the following commutator identity 
for A, B self-adjoint operators 



A 2 B 2 + B 2 A 2 



AB 2 A + AB[A, B] + [B, A\BA + \ ([A, [A, B]]B + B[[B, A], A] 



+ 



A <-> B 



(the second square bracket contains the same expression as the first one with A and B inter- 
changed). After several Schwarz's inequalities, we obtain 

A 2 B 2 + B 2 A 2 > ^AB 2 A + BA 2 B-A 2 -B 2 ^-c[A,B][B,A] (A.51) 
-c[A, [A, B]][[B, A], A}- c[B, [B, A}}[[A, B],B]. 



Using the formula given in Theorem 2.12 (ES-IIIj for the curvature of the covariant derivative 
V = V Q ' a we obtain 



[Dj , D k \ = - V [e . , efc] - d ej d k + d ek dj 
1 3 

-j (e a , Tl{ej, e k )e b )a a a b + ifi(ej, e k ) , 



a,b=l 



where 1Z is the Riemannian curvature, j, k — 1, 2, 3. In short, we can write 

3 

[Dj,D k ] = J2 U ^ D a + W jk , 



a=l 



where Uj k ,Wj k are 2 by 2 matrix valued functions with ||t/jfc||oo < c, ||VZ7, 



< c and 



||W^fc||oo — °b, llVWjfcHoo < cb using (|8.14|) . These estimates guarantee bounds on the double 
commutators as well. 

From these estimates and (jA.51|) it follows that 

D 2 D 2 k + D 2 k D 2 > \ (D ] D 2 k D j + D k D 2 D k ) - cbB 2 - cb 2 . 
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Therefore 



:D 2 ) 2 + (P6) 2 /4 > \ 



E D J + E { D i D l D i + D kD]D k ) 
+ £ (\D* - cbD 2 ) + (P6) 2 /4 - cb 2 . 



(A.52) 



The second line is bigger than (Pb) 2 /8 if e is sufficiently small (P = e 5 ). Every term in the 
first line is nonnegative, so we can complete the estimate ()A.49|) using (jA.50|) and (jA.52|) 



DjD k 



Ab 2 



[Y) 2 ) 2 + (P6) 2 /4 



D k Dj < DjDf. 



8b 2 



D k D 2 D k + (Pby/8 



D k Dj < 8b 2 . 



This completes the proof of ()11.17j) . 

The proof of and are straightforward from (jll-lljl . flll.1511 . (jll.ieij) and 

(110. 14j) . Finally pi.2()j) is proven in the same way as but now directly on the space 

L 2 (dO®C 2 . □ 



A. 8 Proof of Lemma 111 .3t Estimates on the resolvent with a con- 
stant field 

The proof of (jll.21}) - (jll.28|) may be done by straightforward explicit calculations since the 
magnetic Schrodinger operator with a constant field is exactly solvable. We show below how 
to obtain these estimates in a reasonably short way. 

We work in the £ coordinate system, and use that = 0. Because of translation 

invariance in the third direction, we can assume £(u) = 0, so in the £ coordinates we need to 
estimate the operator kernels at (0,0). 

We recall the decompositions (jll.8j) - (jll.l0j) and let 

-a := n 2 + n 2 = {-id, - 16) 2 + H<9 2 + f£i) 2 

be the two dimensional magnetic Laplacian that commutes with n 3 . For simplicity, we denoted 
dj := d%. . By Lichnerowicz' formula ()9.8|) . V 2 = — A + II 2 + a 3 b and recall that n 3 = — id 3 . 
The key idea is that the heat kernel of A has a closed form (see, e.g., Chapter 15 in S79 ) 



47T smh(ot) 



bcothjbt) 2 lb 

^ \Kl — c,±J - jlUi-Wa 



(A.53) 
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Then the resolvent can be expressed as 

1 



V 2 + P 

We define the following norm on R 3 



e -t(P+*H) e tA e td$ dt 



(A.54) 



{btl+Ptl) 1 ' 2 . 
Lemma A. 10. Let P < cb, then the following bounds hold 

1 , e -cK-C\t 

1 - - -7-1/2 e 



V 2 + P 

n 3 



v 2 + p 
ff 



(to 
(to 



< C6P- 1 / 2 



< c6 



e -c|f-CI 



£> 2 + P 



£> 2 + P 



' (to 

(to 



< ce -«-CI 



-1/2 e 



u - CI 2 

1 



J = 1,2 



+ 6 



(A.55) 
(A.56) 
(A.57) 
(A.58) 



where \\ ■ \\ refer to the 2 by 2 matrix norm of the operator kernel as a function from R 3 x R 3 
into the set of 2 by 2 matrices. 

Remark. If b < cP, then the same estimates (|A.55J) - (|A.58J) hold with b replaced by P 
everywhere, including in the definition of ||| • |||. 

Proof. From ()A.53|) and (| A. 54(1 we estimate 

1 



V 2 + P 



(to 



< cb 



oo tP+tb 



Vtsmh(bt) GXP ^ U 



6tcoth(k)(a-a) 2 + (6-C3 



dt 



We split the integration into two regimes: bt < 1 and bt > 1. In the first regime we use 
btcoth(bt) > 1, sin(6t) > bt. In the second regime we estimate coth(fri) > 1 and sinh(6t) > 



\e bt . We obtain 



V 2 + P 



(tO 



< c 


IL 







'1/6 e -tP 



t 3/2 



exp 



dt 



+bcx V ( - i (a-c±) 2 



oo e -tP 



, fV 2 
1/6 1 ' 



exp - 



(fr - C 3 ) 2 

4t 



dt 



74 



after extending both integrations over (0, oo) and using the resolvent kernels of the one and 
three dimensional free Laplacians. The proofs of ()A.56|) - (jA.57|) are similar and left to the 
reader. 

For the proof of (|A.58[) . explicit calculation and trivial estimates yield 



V 



V 2 + P 



< c 



oo e -KP-^ t {cath{bt) - 1)1^ - 



t 3 / 2 sinh(fa) 



bcoth(bt) , 
xexp( -^|£j 



< c 



1/6 P -t(P+b) 



< c 



t 5 /2 

+cb e- cb ^~^ 2 

1/6 e -t(P+b) 



oo e - tP 



-ai 

exp 

la - o 



(6 - c 3 ; 

(£ - C) 2 



o 



f 2 



1/6 

exp 



t 3/2 



4f 
exp 



df 



- Cs) 2 
4f 



+c6 



1/2 C 6| €x -Cx| 2 



oo e -tP 



1/6 



t 3/2 



f 

exp 



df 



(6 - c 3 

At 



dt 



< c e 



-v / E+i 5 IC-CI 



< c e 



t^exp 

+c6 e " cf>l5± " ?±|2 " cv/p|6_?31 

1 



\t-C\ l 



□ 



dt 



dt 



With the estimates of Lemma IA. 101 at hand, the proof of Lemma 111. HI is straightforward. 
For example, the proof of (jll.21|) is as follows 



tr 



(V 2 + P) 



-(0,0) 



R3 



' (0,0 % Mi<<tfp- 1 



V 2 + P 



R' 



-df < cbP- 3/2 



after a change of variables. The other inequalities are proved similarly. □ 
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